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^ I The low temperature expansion of the free energy in a Casimir effect setup 

O^- is considered in detail. The starting point is the Lifshitz formula in Matsubara 
representation and the basic method is its reformulation using the Abel-Plana 

I formula making full use of the analytic properties. This provides a unified de- 

■ scription of specific models. We re-derive the known results and, in a number 



includes the derivations and should provide a self contained representation. 
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1 Introduction 

The Lifshitz formula is the basic tool for the calculation of van der Waals and Casimir 
forces between two material half spaces. It emerged in 1956 for the description of the 
electromagnetic dispersion forces. Together with Casimir's approach of zero-point or 
vacuum fluctuations these are two sides of one coin. In the language of quantum field 
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first interface 



second interface 



Figure 1: The configuration of two interfaces separated by a gap of width a. 

theory these are one-loop corrections to a classical background which may be given by 
boundary conditions or by classical fields as well. 

For the configuration of two parallel interfaces with a gap of widths a between them 
(see Fig. 1), the Lifshitz formula 

^-'■"^T±' f^j:in{l-r,r,e-"") (1.1) 

l^Q J \ I TE,TM 

provides the separation dependent part of the free energy of the electromagnetic field 
at temperature T in terms of the refiection coefficients rj (i=l,2) of the two interfaces. 
In (1.1), the prime on the sum indicates that the (Z = 0)-contribution must be taken 
with a factor 1/2. Initially this formula was written for dielectric half spaces with 
permittivity e behind the interfaces with the well known reflection coefficients 

7] — K eri — X 

r-T-E = — : , = : , (1-2) 

rj + >c erj + K 

which must be inserted for according to the polarization. In (1.1), the integration is 
over the wave numbers k = {^1,^2} in the directions in parallel to the interfaces, the 
summation is over the Matsubara frequencies 

= 2Tik^Tl {I integer) (1.3) 

and the notations 



V = \lil/c' + k\ X = ^ei^/c^ + k\ (1.4) 

[k = |k|) are used. Here, and in (1.1), the permittivity is allowed to be frequency 
dependent, e = e(w), and must be taken at imaginary frequency, u = 
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The Lifshitz formula (1.1) is basic for the calculation of the force F acting between 
the interfaces (actually the pressure since (1.1) is the free energy per unit surface), 



At present, this force can be compared with results from force measurements on a 
high level of precision, see [1] or the recent review [2]. Also, the Lifshitz formula was 
generalized to arbitrary, non-flat geometry of the interfaces in terms of the scattering 
approach, see, for example. Chapter 10 in [1]. 

Along with the great success, it must be mentioned that a decade ago a problem 
appeared the Lifshitz formula has when including dissipation. For instance, when 
inserting the permittivity of the Drude model, 

e^^{uj)^l-—^—, (1.6) 

with a temperature dependent dissipation parameter 7(T), decreasing for T — > 0, the 
free energy (1.5) violates the third law of thermodynamics (Nernst's heat theorem) [3]. 
This means, the entropy 

S = -P, (1.7) 

has a non- vanishing hmit for T — > 0. This problem is still unsolved. Another, closely 
related problem appears if decreasing the dissipation parameter 7 at fixed temperature. 
The naive expectation would be to re-obtain in the limit the free energy of the plasma 
model. While the permittivity (1.6) turns into that of the plasma model, the free 
energy has an additional contribution for 7 — )■ 0, i.e., it is not perturbative in 7. 
This is counter intuitive from the physical point of view since a small dissipation 
should have a small effect on the dispersion force. Also, there is growing evidence 
for a disagreement between the predictions following from the Lifshitz formula with 
dissipation and measurements [2, 4, 5, 6, 7]. 

At finite temperature, the vacuum fluctuations appear together with the thermal 
fluctuations of the electromagnetic held. For the latter, one deflnes a characteristic 
temperature, 

which, at room temperature, corresponds to a separation a ~ 4/xm. The free energy 
depends typically on a dimensionless combination, 

2ak^_T_ 
he "Teff- ^ ' 

It is, for most measurements, a small parameter making the low temperature regime 
relevant. 
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At low temperature, the Matsubara sum in (1.5) becomes inefficient since liigfi I give 
significant contributions. Instead, one uses tfie Abel-Plana formula. For a sequence fi 
(/ integer) of numbers, the representation 

°° / r'^ rii 

J2fi= dlf{l)+ ^— ^^(/(^/)-/(-^/)) (1.10) 
1=0 io e -i 

holds, where /(/) is the analytic continuation to the complex /-plane with /(/) = /; for 
/ integer. Thereby it is assumed that the initial sum converges and that /(/) does not 
have poles or branch points for M > 0. In case the derivatives exist, the right hand 
side of eq. (1-10) has an expansion 

/ dim+f'{o)+r{o)+..., (1.11) 

1=0 -^0 

which, if applied to (1.5), gives the expansion for T — > 0. Regrettably, in a number 
of cases, these derivatives do not exist and a more elaborate treatment is needed. It 
should be mentioned that (1.11) was used in the early calculations of the Casimir force 
to perform the limit of vanishing regularization parameter. 

Using the Abel-Plana formula (1.10), the free energy (1.5) can be rewritten in the 
form 

T^Eo + AtT, (1.12) 

where 

Eo = ^,fd^m, (1-13) 

resulting from the first integral in the right hand side of (1.10), is the vacuum energy, 
i.e., the zero temperature contribution, and the second integral, 

fc„ poo -I 

U TE,TM 

with 

$(a;) = i {(p{iuj) - (p{-iuj)) (1-15) 
is the thermal contribution. We introduced the notation 

POO 

<p{0= rffc A; In (1 - nrse-^"") , (1.16) 
Jo 

which will be used throughout this paper. In (1.13) the integration is over imaginary 
frequencies whereas in (1.14) it is over real frequencies and the integrand involves the 
Boltzmann factor i/^e'^/ksT _ ly 

Using representation (1.10) or (1.11), the low temperature expansion was calculated 
in nearly all cases of interest at least in the leading order. In the present paper we give 
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a detailed representation of this expansion in all cases of interest based on eq. (1-14). 
We make full use of the analytic properties of the function (p{C), (1-16), and are able 
to improve a number of known expansions. Using this method we also reconsider the 
derivation of the contributions violating the third law of thermodynamics as well as 
the non-perturbative contribution appearing for vanishing relaxation parameter. 

It is the aim of the present paper to review the low temperature expansion to the free 
energy for the basic models and represent them unified in the framework of eq. (1.14). 
Thereby we consider the asymptotic expansion of the free energy as given by the Lifshitz 
formula, (1.1), for T — )■ 0. This means that we do not consider any corrections which 
are exponentially small and we also do not discuss the applicability of the expansion 
to the one or other situation. On the other side, we try to cover all relevant models 
and try to give a self contained representation which includes all derivations. It should 
enable the reader to follow all calculations nearly without consulting other sources. 
Therefore, the experienced reader may find some places too much going into detail for 
which the author asks for indulgences. A part of the calculations is made machinized, 
using a standard tool. We tried to explain all steps in detail, such that the reader 
should be able to repeat those calculations easily. 

In the next section we collect the basic formulas for the free energy and, in a 
number of subsections, the specific models we are going to consider. These are, after 
shortly recapitulating of the ideal conductor dielectric with fixed permittivity 

e as simplest example for a medium. Next is the plasma model as the simplest model 
describing a metal beyond the ideal conductor. Then we add dissipation by considering 
a Drude model permittivity. The next subsection is devoted to an insulator, followed 
by a subsection for dielectric with dc conductivity as another example for dissipation. 
Finally we consider the hydrodynamic model for graphene. In the third section we 
derive the low temperature expansion and the specific representation we are using. In 
section 4 we derive the low frequency expansions for the specific models. In the fifth 
section we collect the low temperature expansions for all models. Conclusions are given 
in the last section. Some calculations are presented in the appendixes. 

Following the theoretical approach of this paper we use units with h = c = ks = 1- 
Throughout the paper Cr{s) denotes the Riemann zeta function and je is Euler's 
constant. 



2 Basic formulas and models 

In this paper, the basic formula is the Lifshitz formula, mentioned already in the 
introduction. We consider two plane parallel interfaces perpendicular to the 2;-axis, 
with an empty gap of widths a between them. On the interfaces, boundary or matching 
conditions are assumed to be given or the space behind the interfaces is assumed to 
be filled with a homogeneous medium or to be empty in case of graphene. Also we 
assume homogeneity and isotropy in any plane parallel to the interfaces. We denote 
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the Lifshitz formula in the form 



•^=2^EV(6) (2.1) 



47r2 

with 



poo 

(a)= / dk k In {1 - rir2e-'') . (2.2) 

We allow for different reflection coefficients ri and r2 on the two interfaces. One has 
to insert according to the model considered and to perform the summation over the 
polarizations in case of the electromagnetic field. In general, eq. (2.1) is valid for 
any field if inserting the corresponding reflection coefficients in (2.2). This formula 
represents the free energy per unit area of an interface. The summation is over the 
Matsubara frequencies (1.3) and for frequency dependent reflection coefficients ri{uj) 
on has to use their analytic continuation to ri{i^i). 

The Lifshitz formula was originally derived from the fluctuations of the electromag- 
netic field in the gap introducing a random field into the Maxwell equations. At once it 
was mentioned that this field is associated with the 'zero point vibrations' [8] . Indeed, 
writing the vacuum energy (1.13) in the form 



Eo^^l d^^iO, (2.3) 





it is seen that this is just the half sum of the excitation of the electromagnetic field in 
the sense introduced by Casimir in [9] after Wick rotation. After that, the step from 
(2.3) to (2.1) follows simply by applying the Matsubara formalism. 
In (2.2), as compared with (1.16), we dropped a factor 2a, 

^-2av _^ ^-v (2.4) 

This is equivalent to putting a = 1/2 in all formulas. In fact, this is no restriction since 
this exponential is the only place where the separation a enters the Lifshitz formula. 
Formally, this can be achieved by the substitution k — > k/{2a). The dependence on a 
can be restored at any time by dimensional consideration. Since the free energy is a 
density per unit area, its dependence is restored by 

and all dimensional quantities entering J-' must be made dimensionless by the factor 
2a. For instance, one has to substitute the temperature by 

T 2aT (2.6) 

and any frequency, like the plasma frequency, or dissipation parameter, 7 or a, by 

uUp — >■ 2aujp. (2.7) 
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In fact, there are no other parameter to be restored in this paper. At this place we 
also mention how to restore the fundamental constants. For the free energy, in place 
of (2.5), one has to substitute 

(2.8) 

(this is an energy divided by an area) and 

T^^, a;,^^ (2.9) 

he c 

with T to be measured in Kelvin and cUp in 1/s. 

In the configuration of two parallel interfaces, considered in this paper, the polar- 
izations of the electromagnetic field always separate and are commonly chosen as trans- 
verse electric (TE) and transverse magnetic (TM) modes. The corresponding scalar 
amplitudes, <l>(t,x), satisfy the Maxwell equations and, on the interfaces, boundary or 
matching conditions. These amplitudes can be chosen as linear combinations of plane 
waves, 

$(i, x) ~ e-^'^^+'^^Xs), (2.10) 
with a dependence on the coordinate in perpendicular to the interfaces, 

{::::: '^^^o"<::<;r'' '^-"^ 

Here, u is the frequency, k = {ki,k2} are the momenta (or wave vectors) in direc- 
tions parallel to the interfaces, and k^ resp q, are the momenta (or wave vectors) in 
perpendicular direction outside, resp. inside, the gap. The dispersion relations 

eu"^ = k"^ + kl (outside the gap), 
u"^ = k"^ + q^ (inside the gap), 

follow from inserting (2.10) into the Maxwell equations. Here, the permittivity is 
allowed to depend on the frequency, e = e{u) and we use the notation k = |k|. For non 
transparent boundary conditions like Dirichlet or ideal conductor conditions, the waves 
outside the gap do not contribute to the free energy (2.1). For transparent boundary 
conditions like in the hydrodynamic model one has to put e — 1 outside the gap and 
the corresponding momenta are equal, — q. 
After Wick rotation, 

00 = 1^, (2.13) 
one has always an imaginary wave vector 

q^ir] (2.14) 

inside the gap. The wave vector ks outside the gap may remains real or it may become 
imaginary, 

k3 = ix. (2.15) 
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The dispersion relations turn into 



^(^OC^ = — /c^ + (outside the gap), (2 16) 

= —k'^ + rf' (inside the gap). ^ ' ' 



We will use notations (2.14) - (2.16) throughout the paper. 

In the Lifshitz formula, different choices of the integration variables are possible. In 
(2.1), these are k and ^. In this case, one needs to express all other in terms of these. 



= + k^- (2.17) 

Below, the integration over k in (2.2) will be changed for 77, 

/oo 
d77 77 ln(l-rir2e-^). (2.18) 

In that case one has to express 



k = _ ^2 (2.19) 

with the range r] e [^,00). 

As mentioned above, the properties of the interacting bodies enter the Lifshitz 
formula only through the reflection coefficients. For the interface at 0:3 = a, i.e., for 
the right one if looking on Fig. 1, the corresponding mode function is 

*(^3)-| ^g,fe3(x3-a) {X^>a), ^^-^"^ 

where r and t are the reflection and transmission coefficients. These are to be deter- 
mined from the boundary or matching conditions in — a. The reflection coefficient 
r, determined this way, must be inserted for r2 in (2.2). The reflection coefficient ri 
follows, accordingly, from the scattering from the right on the interface in X3 = with 
the choice 

$(x3j - I ^_,^,3 ^ ^ ^,,,3 ^ {I.ll) 

for the mode function. We mention that the reflection coefficients ri and r2 are inde- 
pendent one from another. The Lifshitz formula stays correct for any combinations. 
In case of non-physical choices, like one from the TE and the other from the TM 
polarization, there would be, however, no physical realization for. 

In the remaining part of this section we specify the models which we will consider. 
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2.1 Ideal conductor 



For ideal conducting surfaces, the boundary conditions are Dirichlet for the TE polar- 
ization and Neumann for the TM polarization. The reflection coefficients are 

rxE = -1, rxM = 1, (2.22) 

resulting in equal contributions from the two polarizations to the free energy. In this 
case, the mode functions $(0:3) terminate on the interfaces. Equivalently, one may put 
t = Om (2.20) and (2.21). 

2.2 Fixed permittivity 

For two dielectric half spaces with fixed permittivity e, the reflection coefficients are 

Q-h eq-k3 

in terms of the real wave numbers related by (2.12) with the frequency 00. In terms of 
the imaginary wave numbers we note 

tte = — \ — , ^TM = \ — , (2.24) 

rj + K ST) + X 

and the relation (2.16) applies, including 

X = ^/{e-l)e + V^■ (2.25) 

In this form, the reflection coefficients enter the Lifshitz formula (2.1) through (2.18). 

The reflection coefficients (2.22) of ideal conducting interfaces follow from the above 
with (2.23) or (2.24) in the limit £ — )■ 00. However, this does not imply that the 
free energy behaves the same way. This can be seen already in representation (2.1). 
Consider / = 0, i.e., the lowest contribution to the Matsubara sum. From eq. (1.3) we 
have ^0 = and with (2.25) we note x = 77 in this case and the reflection coefficients 
become 

rTE = 0, rTM = ^^ (^ = 0). (2.26) 

e + 1 

Hence, in the limit e 00, the (/ = 0)-contribution of the TE polarization does not 
deliver any contribution to the free energy while all other contributions deliver the 
corresponding ideal conductor contributions. This behavior was first observed in [10] 
and motivated 'Schwinger's prescription' to take the limit e — )■ before putting I = 0. 
Also, the question on whether this single mode can influence the result much can be 
answered quite easily. In the high temperature limit, the (/ = 0)-contribution delivers 
the leading order contribution and with the vanishing TE contribution half of the 
result is missing. The physics behind this behavior is transparent. A fixed permittivity 
implies a dielectric material keeping its properties at all, including highest frequencies, 
which, of course, does not happen in physics. 
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2.3 The plasma model 



The plasma model appears if one considers the whole space, or a half space behind 
an interface, being filled with a charged fluid (electrons, for example) coupled to the 
electromagnetic field while the half space before the interface is empty. Eliminating 
the dynamical variables of the fluid from the equations of motion (or, in a functional 
integral approach, integrating them out), one comes to the same reflections coefficients 
(2.23) or (2.24) as above where one has to insert the frequency dependent permittivity 
of the plasma model, 

sP\uj) = 1-% (2.27) 

CO 

Here cUp is the so-called plasma frequency and from (2.19) or (2.25) we note 

K= ^ujI + t]-^. (2.28) 

The spectrum, i.e., the mode content, of this model is well known. A recent discussion 
in the context of vacuum energy was given in [11] and here we mention only that 
it is the same as for fixed permittivity with, in addition, a surface mode in the TM 
polarization. This is a mode with real frequency ui, propagating on the interface, and 
decaying exponentially on the vacuum side of the interface. 

The plasma model describes some basic properties of the electrons in a metal. 
Typical values of the plasma frequency are approximately equal 8-9 eV. Its inverse is 
the skin depths 

5=—, (2.29) 

which is 22-25 nm. The reflection coefficients for ideal conducting boundary conditions 
can be obtained from (2.27) in the limit oOp — )■ oo. Also the free energy of ideal conduc- 
tors is recovered in this limit due to the sufficiently fast decrease of the permittivity 
for large frequencies. 



2.4 The Drude model 

This model is an extension of the plasma model allowing for dissipation. Physical rea- 
sons may be Ohmic losses or scattering of the electrons on the lattice or on impurities. 
These are accounted for by a phenomenological dissipation parameter 7, entering the 
permittivity, 

e''^{u) = 1 - (2.30) 

of the model. In the formal limit 7 — ?■ one recovers the permittivity ^^'(cj), eq. (2.27), 
of the plasma model. The corresponding free energy does not follow in this limit. 

At the moment it is not clear whether the use of the Drude permittivity in the 
Lifshitz formula gives correct results or not [1, 2]. Since we do not enter this discussion 
in the present paper, we take this model as is and make only a few comments. 
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The permittivity (2.30) is complex. Being inserted into the Maxwell equations 
(2.12), a non vanishing imaginary part of the frequency results. For 7 > 0, which 
one needs to assume, this describes dissipation of energy. This is in accordance with 
the intention of the model describing losses, finally resulting in heat. The model is in 
accordance with causality and the permittivity; e^^{u), eq. (2.30), obeys the Kramers- 
Kronig relation. This model, taken alone, has no unitarity and a mode expansion of 
the free energy in terms of real frequencies is not possible. In line with this, it must 
be mentioned that after Wick rotation, which is possible for 7 > 0, the permittivity 

is real delivering with (2.1) a real free energy. In this way one obtains an easy-to-use 
formula. Derivations of this procedure were given in [12] and, recently discussed, for 
example, in [13]. 



2.5 Insulator described by oscillator model 

The response of insulators to the electromagnetic excitations is, beyond a fixed per- 
mittivity, frequently described by the permittivity 

N 

£--i-(^) = 1 + y f^— (2.32) 

of an iV-oscillator model, where Uj are the oscillator frequencies, gj their strengths and 
7j > their damping parameters. Here one excludes the case ouj = since that would 
rather to be described by a plasma or Drude model. In the low frequency limit, a; — > 0, 
one comes to a constant permittivity, 

N 

e{0) ^l + J2iL^so, (2.33) 

j=i 

as considered in subsection 2.2. 

It must be mentioned that (2.32) includes also dissipation processes like the Drude 
model. However, the free energy calculated from eq. (2.1) is real and for all non- 
vanishing oscillator frequencies uij ^ this model is not known to have problems [1]. 



2.6 The case of dc conductivity 

At non zero temperature, dielectrics posses, as a rule, some conductivity due to dissi- 
pation processes like in the Drude model. These are, in the simplest case, accounted 
for by an additional contribution to the permittivity £™^"^ (a;), eq. (2.32), 

£^^^(0;) - £'°^"^-(cu) + i^, (2.34) 
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where a is the static conductivity resulting in a dc current. Being inserted into the 
Lifshitz formula, eq. (2.1), this model results in a real free energy. This conductivity 
is typically a function of temperature, c(T), vanishing at T — j- 0. This model has 
problems similar to that in the Drude model mentioned in Subsection 2.4 [1]. Below 
we consider for completeness also the case of a fixed cr, although it may be physically 
less interesting. 

2.7 Hydrodynamic model for graphene 

In the hydrodynamic model one assumes a charged fiuid, like in the plasma model, 
but confined to a plane, i.e., being 2-dimensional. Again, eliminating the dynamical 
variables of the fiuid, the Maxwell equations appear obeying on the pane matching 
conditions. These result in refiection coefficients 

In the mode expansion (2.20) and (2.21) one has to put = q since from both sides 
of an interface wc have empty space. Accordingly, from the Maxwell equations (2.12), 
only the second applies. 

The matching condition for the TE mode is equivalent to a scalar field with a 
repulsive delta function potential on the interface obeying the wave equation 

+ 2^p5(^3)^ Hxs) = q'^^xs). (2.36) 

For the TM mode, the corresponding scalar problem can be formulated in terms of a 
(^'-potential. 

The mode content of this model is quite similar to that of the plasma model con- 
sidered in subsection 2.3. For instance, there is, for each interface, a surface mode. 
In the limit of infinite plasma frequency, cUp ^ oo, the refiection coefficients turn into 
that of an ideal conductor, eq. (2.22), and the free energy of this model turns into that 
of ideal conductors. 

This model was first considered in [14]. In [15] and subsequent papers it was used 
to describe the vr-electrons of graphene and Ceo- It provides a quite good description 
of their properties in interaction with electromagnetic fields at large frequencies. For 
small frequencies the Dirac model [16, 17, 18] provides a better description. 

3 Low temperature expansion for the free energy 

We take the Lifshitz formula in Matsubara representation, eq. (2.1), as starting point 
for the low temperature expansion. The convergence of the sum in (2.1) and of the 
integration over k in (2.2) come from the exponential factor 



e 



-2«V^F+^ (3.1) 



13 



(we restored, for a moment, the dependence on the gap's width a) making especially 
the sum over / fast convergent. This picture changes with decreasing temperature T 
since / enters through the Matsubara frequency C,i = 27rT/, eq. (1.3). Obviously, for ^/ 
becoming large, numbers I > 1/T must be accounted for. Thus, for decreasing T, the 
convergence slows down and equation (2.1) becomes, in the limit, unusable. 

A way out can be found if an analytic continuation of (fi{Ci) to non-integer, in 
general complex, ^ can be found. This gives the possibility to define ip{$,) as a function 
in the complex ,^-plane and, using the Cauchy theorem, to represent the Matsubara 
sum in (2.1) as an integral, 

^EV(6) = |¥^(0) + T f rf/^— i-^^(27rT/), (3.2) 
1=0 "^^ 

where the path T encircles the non negative integers, / = 1,2,..., and crosses the 
real axis in I — 5 with < 5 < 1. The next step is a deformation of the integration 
path towards the imaginary axis. For > 0, i.e., on the upper half of the path, one 
substitutes 

with u) e [0, oo) and the exponential in the denominator becomes large for uj — > oo. 
For < 0, i.e., on the lower half of the path, one substitutes 

Since, in this case, the exponential does not grow for large cu, one needs to rewrite it, 

' 1 - (3-5) 



^ g— i27r( ]^ gi27r/ ' 

In the contribution from the first term on the right hand side it is meaningful to change 
the integration variable according to / = and to write down this contribution 
separately. The integration can go along the real axis since there are no poles in this 
contribution. The second term can be joined with the contribution from the upper half 
of the path. In both cases uj runs from zero till infinity. One comes to the representation 



(3.6) 



This is the well known Ahel-Plana formula. 

In moving the integration path V towards the imaginary axis and performing the 
limit (5 — )■ 0, from the origin, i.e., from / = 0, a contribution appeared which just cancels 
the first term in the right hand side of eq. (3.2). We mention that in eq. (3.6) there 
is no pole for a; = due to the compensation in the parentheses. Further we mention, 
that in (3.6) it is assumed that the function ip{^) is continuous in ^ = 0. If this is not 
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the case, one cannot move the path completely to the imaginary axis. However, such 
situation does not appear in the examples considered in this paper. 

A further assumption in deriving eq. (3.6) concerns the function ^p{^). It is assumed 
that it does not have poles or branch points in the half plane > 0. Otherwise, from 
moving the path there would be additional contributions. This property is always 
guaranteed if the modes of the electromagnetic field are subject to an elliptic scattering 
problem. It holds also for the model with dissipation where the corresponding poles 
are all located in the half plane < 0. For vanishing dissipation parameter, these 
approach the imaginary axis from the left and the path must pass them on the right 
side, for instance, by adding an infinitesimal amount, 

(p{±iuj) (p{±iuj + 0), (3.7) 

which is necessary for all models anyway. 

In apphcation to the free energy (2.1), eq. (3.6) defines a split, 

T^Eo + AtJ^, (3.8) 

into the vacuum energy Eq (2.3), 

47r2 J, 



Eo = :^l d^m, (3.9) 



resulting from the first term in the right hand side of (3.6), and, from the second term, 
the temperature dependent part, 

involving the Boltzmann factor l/(e'^/^ — 1) and 

^uj) ^ i {(fiiio) - ip{-iu})) . (3.11) 

The integration variable u has the meaning of the frequency as entering eq. (2.12) 
in the case a mode expansion makes sense. As already mentioned, this is not the 
case for models with dissipation (see the remark at the end of this section). However, 
independently on the interpretation, representation (3.8) with (3.10) and the property 
(3.11) are valid for these too. 

At this place, an important remark on the direction of the contour rotations is in 
order. The rotation (3.4) is the inverse of the usual Wick rotation (2.13), whereas (3.3) 
is the inverse of an Anti-Wick rotation. Since it is customary to write the Abel-Plana 
formula just with the order of terms as in the parentheses in (3.6) with (p{iuj) first, in 
application to the free energy (3.11), the term corresponding to the inverse Anti-Wick 
rotation, goes first. Of course, this does not change anything except for notations. 
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Below we will find it convenient, in a number of occasions, especially after some 
variable substitutions, to use the refiection property. 



^{-iuj) = ip{iuj)\ (3.12) 

this function with reflection coefficients fohowing from a scattering problem has, to 
represent the difference in (3.11) in the form 

i{ip{iuj) — (fi{—iu!)) — —{(f{iuj) — c.c), (3.13) 

where c.c. denotes the complex conjugate of what is in front to be inserted. In doing 
so one has only to pay attention to signs in some places, especially in the permittivity, 
which changes under Wick rotation e{u!) — >■ s{i^), but which enters (f{iuj) after Anti- 
Wick rotation, 

e{i0^e{-u;). (3.14) 

This sign shows up in models with dissipation only. 

In some simple cases it is possible to use the Abel-Plana formula, formally not 
entering the complex plane. Assume the function ip{^) has a Taylor series expansion, 

^(0 = E^^^"^(0)' (3-15) 

n>0 

one gets for $(0;), eq. (3.11), an expansion directly in terms of real quantities. This 
can be inserted into (3.10). Interchanging the orders of integration and summation, 
the integration can be carried out. One obtains 



At J" = At y (-1)^+'2Cr(2A; + 2) (^(2fc+i)(o) 

k>0 

with the Riemann zeta function in even integers. 



(3.17) 



in terms of the Bernoulli numbers -B„. Obviously, this is an expansion for T — )■ 0. For 
most systems, however, the function ip{C,) does not have a Taylor expansion. Typically, 
in the examples considered in this paper, (p{0) exists, but not the derivatives. Never- 
theless, even in the case there is no Taylor expansion, the derivatives, as far as they 
exist, give with eq. (3.16) the lowest contributions to the asymptotic expansion for 
T ^ 0. 

We add a remark on the convergence of the vacuum and the free energies. In gen- 
eral, the vacuum energy, and with it also the free energy, have ultraviolet divergences 
resulting from slow convergence for large frequencies or momenta. In the situation of 
a Casimir effect setup, considered here, these divergences do not depend on the width 
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of the gap and the Casimir force is always finite. The spht (3.8) is, of course, vahd 
beyond the Casimir effect setup. In that case the divergences in the free and in the 
vacuum energies are the same and the thermal part AtJ^, (3.10), does not have any 
divergencies. Its convergence follows, obviously, from the Boltzmann factor, whereas 
before the contour rotation, i.e., in eq. (3.2), this factor is bounded. It is just the 
contour rotation, which, without changing the integral, redistributes for small T the 
main contribution to the integral towards small u. As a result, for T — )■ 0, the integral 
over u is fast convergent and the contributions from > determine the asymptotic 
expansion for low T. This is in opposite to the situation in the initial Matsubara 
representation where large imaginary frequencies were needed for. 

At this place we mention the Poisson re-summation formula which is yet another 
way to redistribute the convergence. In order to use that formula either one has an 
explicit representation, typically an Gaussian exponential, or one needs to make an 
analytic continuation from integer / to, at least, real ones. One obtains, in place of the 
Matsubara sum, another sum, which is fast converging for small T. We do not use this 
approach in the present paper. 

Due to the convergence properties just discussed, cq. (3.10) allows, in a simple 
way, for the low temperature expansion of AtJ^- For this, it is sufficient to assume the 
function $(a;), eq. (3.11), has an asymptotic expansion for a; ^ 0, 



(w) = (^$1 + l»i Ina;) a; + $2<^^ + ^3<^^ + . . . . (3.18) 



Here we allowed for a logarithmic contribution in the first order since it will appear 
below in the Drude model. If the low frequency expansion (3.18) is found, the asymp- 
totic expansion, for T — )■ 0, of the free energy can be easily written down by inserting 
(3.18) into (3.10) and using 



Jo ^^e-/^-l 



du—j^— = r{s + 1)Cr(s + l)T^+' (3.19) 



for the integration over u. For logarithmic contributions one may take the derivative 
of this formula with respect to s. In this way one comes to 



At J" = ^ 



47r2 



;<^i + (CR(-l)+ln(27rT))$i)CR(2)T2 

+ 2Cr(3)$2T3 + 6Cr(4)$3T^ + ...]. (3.20) 



By using the explicit values (3.17) and multiplying out the square bracket, this formula 
can be rewritten in the form 



At J- = ^ (^1 + (Cr(-1) + ln(27rT))4.i) + ^ ^^T' + ^ ^sT^ + .... (3.21) 

It should be mentioned that expansion (3.20) can be derived for any model if $(a;) 
does not depend on T. It starts always at least from T^. For this reason, it cannot 
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Figure 2: The complex plane of the wave number 77, integration pathes shown. 

come in contradiction with thermodynamics. It should be mentioned that expansion 
(3.16) is a special case of (3.20) since it has only odd powers of T and no logarithmic 
contributions. 

With representation (3.20) resp. (3.21) at hand, the 'remaining task' is the calcu- 
lation of the coefficients $i. For this we return to eq. (2.2), 

(^(^) = / dkk\n{l- rirae"") , (3.22) 

assuming the analytic continuation to the complex ^-plane is done. Following from 
eqs. (2.12) and (2.16), the variables w, k and t] are related by 



7] = VP - uj^. (3.23) 

It is possible to change the integration in (3.22) from k to t], which runs from t] = iu 
till infinity parallel to the real axis. In Fig. 2 this integration path is denoted by Fi. 
In deriving the expansion for the various models in the next section, it turns out to 
be convenient to change the integration path for the sum of two, one running from 
7] = iu to 1] = 0, and the other, from rj = along the real axis till infinity. These 
two are shown in Fig. 2 as r2 and T^. Of course, the integral does not change. In 
representation (3.22) with integration over k, this corresponds to a subdivision of the 
integration region into two regions. 



region (a): < k < u, with q = y/u"^ — k^, , . 

region (b): to < k, with 77 = ^k'^ — w^, 

where, at once, that wave numbers are shown which are real in the given region, as 
it follows from eqs. (2.12) and (2.16). As a convention, we will these regions always 
denote by the letters a and h in parentheses. 
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According to (3.24), the integral in (3.22) splits into two, 

ip{iuj) = ip(a){i<^) + V^{b){ii^): (3.25) 
which will be treated separately. With eq. (3.11), this induces a corresponding spht 

$(a;) = $(„)(a;) + $(b)(^), (3.26) 

and the relations 

$(6)(a;) = i(93(fe)(ia;) - v3(6)(-ii^)). (3.27) 

In region (a) we have 

(^(„)(ia;) = / dkk]n{l- rirae"^*) (3.28) 
Jo 

with q shown in (3.24). The reflection coefficients must be expressed in terms of oj and 
k. For instance, for the wave number wc note 



ks = Veu^ - k^, (3.29) 

which may be both, imaginary or real in dependence on the model. Changing for the 
integration variable q, eq. (3.28) can be written in the form 

<f^a) (ioo) = / dqq\n{l- nr^e-"^) , (3.30) 
Jo 

where one needs to express everything in terms of q and u, for instance k — ^^u^ — q'^ 
and ks — a/(F^^T)cu2~+~^. 
In region (b) we have 



f oo 



/■oo 

/ dkkhi{l — rir2e~^') 

J O) 

/■oo 

= / dry 77 In (l — rir2e~'') . (3.31) 
Jo 

In the second hne we changed the integration variable for 77 using (3.24). This inte- 
gration corresponds to the path in Fig. 2. The reflection coefficients entering the 
second line must be expressed in terms of u) and 77. For instance, from (2.16) we note 



x= V-(£- l)w2 + ^2_ (^3 32) 

We remind that we keep, in all calculation, the relations x = ikz and 77 = iq. 
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4 The low frequency expansion for specific models 



In this section wc obtain the low frequency expansions (3.18) of the function ^{u) for 
various models. This section comprises the main technical part of the paper. Some 
calculations are banned to the appendixes. As mentioned in the Introduction, we 
use the simplest form of notations, especially we drop the factor 2a everywhere as 
announced. 



4.1 Ideal conductor 

This is the simplest model and well known. We consider it for completeness. At once 
we illustrate the technique used, especially the division of the integration in (p{iuj) into 
two regions. Also it allows for an easy checking of the overall factors. 

The reflection coefficients for ideal conductors are given by eqs. (2.22) and their 
product, entering (p{iuj), is rir2 = 1 for both polarizations. The contribution from 
region (a) can be written in the form 

iPia){ico)= / dqq hi {l-e-"'), (4.1) 

where we changed the integration variable for q using (3.24). The logarithm can be 
written in the form 

ln(l-e-^^) = -| + ^|+ln(2sin|). (4.2) 

Now, as long as a; < tt, the sine does not change sign and the and the logarithm in the 
right side stays real. Now we calculate according to (3.27) the imaginary part, 

^(a){u)) = i(v?(a)(«w) - C.C.), 
/•oo 

= / dqq{q-7r). (4.3) 

^0 



The remaining integration is trivial. 



In region (b) we note 

/•oo 

ipQ,){iu:)^ j dr]r]\n{l-e~'') . (4.5) 
Jo 

This expression is completely real. Hence it does not contribute, $(b)(a') = 0, and we 
get from (3.26) and (4.4) 
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We mention that this formula not only provides the asymptotic expansion for a; — > 0, 
it is exact for < tt. 

From eq. (4.6), in the context of eq. (3.18), we have non vanishing coefficients 

TT 1 

*2 = ^3 = g- (4.7) 

Inserted into eq. (3.21), these deliver the expansion for ideally conducting interfaces, 

Aj.r" ^ -^^T' + —T^ + . . . . (4.8) 
27r 90 ^ ' 

Here we included a factor of 2 to account for the two polarizations of the electromag- 
netic field. The dots represent exponentially decreasing contributions we do not care 
of in the present paper. 

It should be mentioned that the method used in this paper is only one out of quite 
a number of equivalent ones applicable for this simple model. More details can be 
found, for example, in [1], chap. 7.4. 



4.2 Fixed permittivity 

Dielectrics with fixed permittivity e represent the simplest model for an insulator. 
While the ideal conductor considered in the preceding subsection is in the aim of 
Casimir's original idea, an insulator is rather in the spirit of Lifshitz's approach. Here 
both are treated within the same formalism. Also, the model with a fixed e may serve 
as a good approximation for more complicated permittivities at low frequency. 

The reflection coefficients are given by eqs. (2.23), or by (2.24), which will be used 
in regions 1 and 2, accordingly. As already mentioned, the limit £ — )■ oo does not 
turn the free energy into that of ideal conducting interfaces. Hence, a situation with 
one interface ideal conducting, the other with finite permittivity e behind, is different 
from a situation with both interfaces having finite permittivities £i and £2 behind and 
cannot be obtained by any limiting process. For this reason we consider 4 cases as 
shown in Table 1 and denote the case as an index in parentheses, $(A;)(a') (k=l,...,4). 



case 


first interface 


second interface 


polarization 




ri 


r2 




1 


-1 


tte 


TE 


2 


^TE 


vte 


TE 


3 


1 


ttm 


TM 


4 




fTM 


TM 



Table 1: Notations for the four cases considered. 



As discussed in Sect. 3 we will perform the calculations separately in regions (a) 
and (b). We add this information to the index such that 

^(fe,n)H (^=l,-,4, n^a,b) (4.9) 
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denotes the contribution from region n to case k. We use his notation also for the 
functions ip{iijj) and in the relations 

*(fc,n)('^) = i {V{k,n){i'^) - c.c.) . (4.10) 
For a given case, the contributions from both regions must be added, 

In the final result for the electromagnetic field, cases 1 and 3 or cases 2 and 4 must be 
added. 

4.2.1 Region (a) 

Here we use the reflection coefficients as given by cq. (2.23). In eq. (3.28) we change, 
for convenience, the integration variable for q = \/ixP- — k'^. In this case, in eq. (2.23), 
one has to use k^ — -^(e — l)c<;^ + and we get 

^{k,a){iuj) = dq q In {1- rir^e'''^) , (4.12) 

^0 

where for ri and r2 one has to insert according to Tabic 1. 

In this expression, a direct expansion of the integrand in powers of u with subse- 
quent integration over q delivers the expansion of ip(^k,a){i'^) and, by means of eq. (4.10), 
that of $(fc,a)('^)- In the latter only odd powers of uj remain. Defining expansion coef- 
ficients $(fc,a),i(a;) in parallel to eq. (3.18), these can be calculated easily by machine. 
The non vanishing coefficients are shown up to the order i — 6, which corresponds to 
an order in the expansion (3.21), 



$(l,a),3 



li-e'/' + {e-lf' + l), (4.13) 



%a),3 = ^((Vi+l)(Vi-V^)-2), 

%a),5 = ^|^(3yi(x/i-V^)-l), 

%a),3 = ^^^((l-2(v/i- (Vi+l)v^-l), 



^(3,a),5 



: ^^-^ (Vi (48£3 - 72£2 + 10£ + 15) - 8(£ - If" (6£2 - l)) , 
180 

%a),3 = {esj^e - 1 {2e'l' + As + ?>^e + 2) - {^e + if [2e' - l)^ , 

*(4,a),5 = fe^(s^/^(2.-l)(24s^-7s-3)+x/^(-48.^ + 8s^-2s^-3)). 



22 



4.2.2 Region (b) 

In this region we use eq. (3.31) and insert the reflection coefficients as given by eq. 
(2.24) with X = -^/rf — {e — l)uj'^, which foUows from eq. (2.19). Next we observe that 
in 

poo 

k,b){ioj) = dr]r]ln{l- rirae"'') (4.14) 
Jo 



the part of the integration with rj > y/e — loo dehvers a real contribution and does not 
contribute in (4.10). Thus we restrict the integration region, 



*(fc,fe) (^) = / dri rii {In {1- nrae"") - c.c.) . (4. 15) 
Jo 

In this expression x is imaginary, x = iks with — {s — l)u;^ — rf. Now we make 
the substitution 77 — >■ y/e — la; 77 and get 



$(,,,)(a;) = (£-l)a;^*(^^^ill^) (4.16) 



with 



= ! dr]r]i{hi{l- nrse-^^") - c.c.) (4.17) 
Jo 



and for the reflection coefficients, in dependence on the polarization, one has here to 
insert 

™ ~ , • n — r ^™ ~ — , • y^-^^i 

For the expansion of the function $(fc 2)(i^); eq. (4.28), for a; ^ 0, we need to know 
the expansion of the function ^'(/x) for /j, ^ 0. 

It is not possible to expand the integrand in (4.17) since a subsequent integration 
would not converge for 77 —> 1 since for the reflection coefficients (4.18) the relations 

limrTE = l, limrxM = 1 (4.19) 

hold, which result in singular terms in an expansion of the logarithm in (4.17). We 
proceed by factorizing the logarithm in (4.17) and representing ^(//) as a sum, 

'"i'"2, A*) + ^_b('"i'"2, A*-) (4.20) 



of two functions. 



*a(^, A*) = / (iry 77 i (In (1 + re '^'') - c.c.) , 
Jo 

^^(r,/x) = / d7;7;i(ln(l-re-^'') -c.c.) . (4.21) 
Jo 
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The functions $(fc,6)(a;), defined in (4.9), can be reobtained from from these as follows, 



easel: $(i^(,)(a') = {e — l)uj^ ^a(1, — 1<^), 

case 2: $(2,6)(u;) = {e - 1)uj^^a{1,IV^^uj) + ^b{1,IV^^uj)) , 

(4.22) 

case 3: $(3,6)(t^) = (e - l)w^ ^b(£, - Iw), 

case 4: $(4,6)(i^) = (e - l)a;2 (^^(e, - Ic^) + *b(£, - 1^^)) , 

where we used rTM|e=i = '"te, which holds for the coefficients (4.18). 

For the function ^^(//), the expansion can be obtained easily by expanding the 
integrand with subsequent integration. The singularities resulting from the expansion 
of the logarithm appear for 77 = and are compensated by factors from the expansion 
of the exponential. We denote this expansion in the form 

^A{pi)^J2'^A,i{e)pi\ (4.23) 

i>0 

The coefficients up to order 4 are shown in Table 2. 

For the function ^'b(/^), which carries the above mentioned singularities if expand- 
ing the integrand, we define an auxiliary function 



e — i^\ — if 



*aux(/i) = C dr]y]i\\Y^\\- - ^^^^ e-^ ) - c.c. ) . (4.24) 



In this function, the integration can be carried out explicitly with subsequent expansion 
in powers of We use this function to represent ^b(//) in the form 

*S(/X) = *sub(/^) + *aux(/^) (4.25) 

with 



*sub()")=/" rfr/r/i (in (l-re-"^) -In (1- ^ — ^% -^|-c.c.|. (4.26) 
h \ \ e^%^\-ri^ ) J 

In this function, the integrand can be expanded up to the order of //^ with convergent 
subsequent integration over r]. In this way we obtain the expansion 

4 

*B(/i) = ^^B,i(£)//^ + .... (4.27) 



1=0 



All these operations can be carried out by machine. The coefficients up to order 4 are 
shown in Table 2. 
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i 





1 


2 


3 


4 


WA,i{£) 




1 










2(£ + l) 


3e 


180£3 




ne 


2s 








2{e+l) 


3 


8 


45 


48 



Table 2: The coefficients WA,i{s) and WB,i{s) appearing in eqs. (4.23) and (4.27). 



Now, using tfie expansions (4.23) and (4.27), we can calculate from (4.22) the 
expansion 

6 

^ikM^) ^J2^ik,b),i^' + 0{uj'), (4.28) 

i=0 

where 

= WA,i{e) + WB,i{s). (4.29) 
For the four cases defined in Table 1, these expansions read 

= lie - - \{e - 1)3/V + l(s - Ifl-u^ + O (^^) , (4.30) 



*(2,6) 



(3,6) 



lie - ifW - ^ie - 1) V + l(s - 1)VV + ^(s - 1) V + O (c^) , 

_7r (£-1)^ 2 2 _ / _ TT _ 

2 (e + 1) 3^ ^ 8^ ^ 

+^ie - l)'/'e i6e' - l) co' + ^(e - 1) W + O {cu') , 



7r7i{e-lf 2 £-l'/'2£2-l vr , ^,^2 4 
$,4^^ = } /-^j^j^l 1 !^ (e-iye^oj^ 



ie - 1)5/2 ^4g^6 _ 8^4 ^ + 3) TT , ,,3 2 6 ^ / 7^ 

In this way, the calculation of the contributions from region 2 to the low temperature 
expansion is finished. 



4.2.3 The low frequency expansion 

Here we collect the contributions from the two regions calculated above. According 
to Eq. (4.11) one has to add them. It is seen that from region (a), eq. (4.13), only 
odd powers of uu come, whereas from region (b), eq. (4.30), even powers come in too. 
Together, the expansions coefficients for the four cases in Table (1) are 

4(„ = - + i (1 ^3 +^(,_ 1)^3/2^5 + (^7) _ (4 3^ 
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{nie - l)e) 



^(3) = — ^ 



2{e + l) 3^ ^ 8 

- 1V3/2 

4 



^ (2e(12£(2£ - 3) + 5) + 15)^^ + -{e - 1) W + O {co') , 



*(4) = -^^^|^^^ + ^(v^(.(2.-3)-l) + 2)a;^-|(e-l)W 

1440 768^ ^ ^ ^ 

For the electromagnetic case we have to add the polarizations. In case of two dielectric 
half spaces we have to add cases 2 and 4 from (4.31), 

+^(e-lf(e' + l)co' + 0{co'). (4.32) 

For one interface ideal conducting in front of a dielectric half space we have to add 
cases 1 and 3, 

^id.cond.,£M = + $(3)(a;), (4.33) 



^~ + ^ (^5/2 _ 2^3/2 + 1) ^3 _ _ ;l)2^2^4 



4(£ + 1) 3 ^ ^ 8 



, ,5(12£(2£ - 3) + 5) + 10)0;^ + -{e - 1) W + O {u') . 

Inserted into eq. (3.21) these coefficients give the expansion of the free energy up 

to T^, see Section 5. As aheady mentioned, the hmit e — )■ 00, does not turn these 
coefficients into that if of the ideal conductor. Also, the two cases (4.33) and (4.32) 
are independent from one another. 

4.3 The plasma model 

The plasma model is described by the reflection coefficients (2.23) or, equivalently, 
(2.24) and by the permittivity (2.27), 

e^\iu) = 1-4- (4-34) 
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Since in the limit of infinite plasma frequency, 0;^ — > oo, the free energy turns into that 
of the ideal conducting interfaces, there is no need here to introduce the cases used in 
the preceding subsection. Instead, we use on each interface its own plasma frequency, 
denoted by cupi and U!p2. The case of one ideal conducting interface can be restored 
afterwards by sending one of these frequencies to infinity. 

For the calculations, we divide the contributions according to regions (a) and (b), 
defined in Section 3, Eq. (3.24). 

4.3.1 Region (a) 

In this region we have a real q and we use representation (3.30), 

^(a){i(^)= / rfgg ln(l -rirse-*") , (4.35) 

where we have from (2.12) with the permittivity (2.27) an imaginary = in with a 
real k = 1^0;^ — g^. The reflection coeflicients are given by eqs. (2.23), where one has 
to insert for the corresponding interface, which we indicate by an additional index i 
(i=l,2), 

rXEi = ——, TTMi = 7 -T , (4.36) 

with Xj = Y^cUpj — For sufficiently small cu, uj < cUpi, these reflection coefficients are 

pure phase factors, i.e., their modula are equal to unity. Thus we can rewrite them in 
the form 

TTEi = e-^('^+2"-'), rxMi = e-^2"-S (4.37) 

with 

q ^^J^pi-Q^ 
ctEi = arctan — , , = arctan kt-^- (4.38) 



pi 

This allows to cast (4.35) in the form 



(„)(H= r dqq In (4.39) 
Jo 



where one has to insert 

ipE — 2q;ei + 2q!e2 + ? (TE polarization), 
■0M = 2q;mi + 2q;m2 + ? (TM polarization). 



(4.40) 



Like in the case of ideal conducting interfaces, eq. (4.2), this simple structure allows 
to rewrite the logarithm as 

In (1 - e-''^) = -i^ + i- + In (2 sin ^"j . (4.41) 

2 2 y ^ / 
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As before, the remaining log does not contribute since its argument does not change 
sign for sufficiently small ui and we get with (3.11) 



$(a),TE(tc') = / dqq{i)^-T:) (4.42) 

and 

^ia),TM{oo)^ dqq{ipM-T^)- (4.43) 
Jo 

In these expressions, the integrations can be carried out explicitly and we come to 

^{a),TM{oo) = -^a;' + ^a;=^ + 2a;Ji/iM ( — ) +2cc;2 /iM ( — ) , (4.44) 

I 6 \^pl/ X^p'iJ 



where 



h^{z) — I dqq arctan 
Jo 



j^r^^ 2^2 _ I 

\ arctan z (4.45) 



and 



huiz) = / dqq aictan — J- — 

Jo ~ ^ ) 



H ^ ^ arcsm z (4.46) 



result from the integrations in (4.42) and (4.43). 
The expansion of these functions. 



M^) = J2h^,^''-l'' + ^^' + ^'' + Oiz% (4.47) 

i>3 

hu{z) = J2hM,sZ^-z^-'^z' + lz'-lz'+^^z' + 0{z% 

i>3 

can be used in (4.44) to obtain the expansions 

$(a),x(u;) = -7ra;2 + + i ("^ + /.^ 3^) o;^ + ^ ( ^ + A_ j h^ .^i^ 

(4.48) 

where we used the notation X for TE and TM not to write down nearly the same 
formula twice. From this representation it is seen that the ideal conductor result (4.6) 
can be reobtained in the limit of both plasma frequencies becoming large. 
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4.3.2 Region (b) 

In this region we have a real 77 and, from (2.25) with (4.34) inserted, and a real x = 
+ yf"- The reflection coefficients (2.24) are, explicitly written, 

rTE = rTM = j ^ ; ■ (4-49) 

In the case, the two interfaces have different plasma frequencies, wc have to insert Upi 
and ujp2 into the corresponding reflection coefficients {i — 1, 2), eq. (2.24)and in 



X,; 



= (i = l,2). (4.50) 

These coefficients are real. Thus the function (3.31), 

/■oo 

</7(b)(ia;) = / dr) r] In (1- rirae"") , (4.51) 
Jo 

seems to be real too. If this would be the case, it would not contribute to 

$(6)(a;) = i(^(6)(ia;) - c.c). (4.52) 

However, it happens that the argument of the logarithm in (4.52) changes sign not 
only for a finite u as in the case of an ideal conductor, but also for arbitrarily small cu. 

Physically this can be understood from the spectrum. We consider real u with 
imaginary both momenta, q = if] inside and = ix outside the gap. The plasma 
model is known to have, for the TM polarization, such excitations, the surface modes 
(or surface plasmons). Further, the argument of the logarithm in (4.51) is proportional 
to the transmission coefficient of the system of two interfaces, thus it has zeros at the 
wave numbers corresponding to these excitations at a given frequency u. Moreover, 
the argument of the logarithm has poles corresponding to the surface plasmons on 
the interfaces taken individually, where their reflection coefficients (i=l,2) become 
infinite. 

Since there arc no surface plasmons for the TE polarization, wc can be sure that 
the corresponding logarithm in (4.51) stays real and that his polarization does not 
contribute to (4.52). For the TM polarization we rewrite (4.51) using (2.24), 

POO 

ifip) (iu) = / dr)r) [-In {eiT) + xi) - In (£2^7 + X2) (4.53) 
Jo 

+ In ((£17; + xi) {627] + X2) - {siT] - xi) (£2^ - X2) e"'')] . 

Here, the arguments of the logarithms have no poles but only zeros. We denote the 
wave number 7] for the single plasmons by 77singie,i (^ — 1)2). These follow from the 
poles of Ti, i.e. these are solutions of 

EiTj + Xi^O (i = l,2). (4.54) 
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Using (4.34) and (4.50), these equations can be solved for 77 resulting in 



^2 



?7single,i = / (« = 1, 2). (4.55) 



u^l, - 2a; 



The wave numbers of the plasmons in the system of two interfaces we denote by 77symm. 
and r/antisymm.- Thcsc notations account for the properties of the corresponding wave 
functions to have a symmetry in case of equal plasma frequencies. These wave numbers 
are solutions of the equation 

1 -nrse"'' = 0, (4.56) 

or, equivalently, of equating the argument of the third logarithm in (4.53) to zero. 

It is known (see, for example Fig. 2 in [11]) that r^antisymm. is not real for small u). 
Therefore it does not contribute in region (b). The solution 77symm. of eq. (4.56) exists 
for arbitrarily small ui and it has an expansion 

?7symm. = ^ % (4.57) 

whose coefficients can be easily calculated by inserting (4.57) into (4.56) and expanding 
in powers of oo. One obtains 



?7symm. = \/2 ( — + — ] 



6 6 6 1 

+ —= I —r + + ^ + ^ + 



6-\/2 \^pl (^pl^p2 l^pl l^pl <^pl<^p2 <^pl<^p2 
'^p2 ^p2 ^p2/ \i^pl+^p2 

In case of equal frequencies, this expression simplifies, 

?7symm.| = 2-^ + ^ ^ ^ ^ 1 \ 



f I 3 7 1 11 \ / w V ^ / 



The analytic continuation in <^(b)(ia;), eq. (4.51), is to be taken starting from large 
77 where the expression is real. When the argument of a logarithm, for decreasing 77, 
passes zero, the logarithm acquires an addendum of in. Thus we have in (4.52) 

/■'7single,l /■»7single,2 /■»?Bymm. 

$(b)(a;) = -27r - / drjT)- dr]r]+ drjr) , 

I Jo Jo Jo 

— ~^ (^symm. ~ ^single,! ~ 'ysingle,!) ■ (4.60) 
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Inserting (4.55) and (4.58), after re-expansion, we get 

$(;.)(a;) = -27r( — + — ) 0;' + ^ (2 (^ + ^ ) +^ 



,^pl ^pl^p2 ^^pl ^pl^p2 ^p2 ^p2, 



which is, for equal plasma frequencies, 

I = -47r— + — 1 ^ ^ + O . 4.61 

It is seen that in the limit of large plasma frequencies, there is no contribution from 
region (b). 

4.3.3 The low frequency expansion 

Adding the contributions from both regions we get the low temperature expansion for 
the plasma model. For the TE polarization, only region (a) contributes to $(0;) and 
we have only to rewrite (4.48) adding the index 'TE', 

For the TM polarization we have to add (4.48) and (4.61), 

V2 yupi u:p2j J V3 V^^pi Wp2/y 

+ (l _ ^^p2 _ 6a;p2 6 ecupi _ Qupi ^ Upi 

3 \ cjpi ujpi ujp2 ^p2 2a;p2 

+^.^)^^in^J_)^^ + 0(.^), (4.63) 
which also simplifies, 

+— 1 -)— + -— + 0{uj'), 4.64 

in case of equal plasma frequencies. 

Finally, for the complete electromagnetic case, we have to add the polarizations, 
eqs. (4.62) and (4.63), and get 

V \i^pi ^p2j J V3 3 \Upi ujp2j J 
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3 V <^i!i w^i (x'^i w^o <x'p2 2a;p2 
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'^-;jr]z^ + f{dr + ^]-' + 0(.'-). (4.65) 



The coefficients of this expansion, identified according to eq. (3.18), and being inserted 
into eq. (3.21), give the low temperature expansion for the plasma model. 
As special case we note the corresponding formula, 

$ed(c^) = -TT ( 1 + — 1 cj^ + - I 1 + — 1 o;^ 

+— 1 2)— + T— + ^ ^ ' 4-66 

for equal plasma frequencies and 

^EDioj),, = -TT ( 1 + — I CU^ + - I 1 + — 1 

'^'-'V /Wp 1=00, 2=<*'p \ ,./ o\ ../ 



y CUp J 3 \ UJp / 

TT / 12 12\ 00^ 12 cj^ ^, g, 
6 V UJp ujU ujI 5 a;3 



for one interface ideally conducting. 



4.4 The Drude model 

The reflection coefficients for the Drudc model are given by eq. (2.23) or by (2.24) 
with the permittivity given by eq. (2.30). As already mentioned, the limit 7 — )■ 
of vanishing dissipation parameter does not reproduce the free energy of the plasma 
model and there is a problem with thermodynamics. In the following subsections we 
consider first the case T — >• for fixed 7 and, afterwards, the case of vanishing 7. 



4.4.1 Fixed dissipation parameter 

In this subsection we consider a fixed dissipation parameter 7. In that case the low tem- 
perature expansion of the free energy is given by eq. (3.21) with coefficients following 
from the expansion (3.18). 

As already mentioned, the limit 7—^0 does not reproduce the free energy of 
the plasma model. Thus, even for fixed dissipation parameter, the low temperature 
expansion of the free energy in the Drude model is quite different from that of the 
plasma model and the ideal conductor results, on one or on both interfaces, and it 
cannot be obtained as some limiting case. For this reason, we consider the four cases 
introduced in Section 4.2. in Table (1). Also we restrict ourselves to the order u, i.e., 
to the order T^, which is the leading order for this model. On the one side, higher 
orders are technically elaborate, on the other side, these are hardly of any use. 
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For the calculations we use the scheme of two regions introduced in Section 3, eq. 
(3.24). In region (a) we use representation (3.30) with q as integration variable, 

' dq q In {1 - nr2e-"') . (4.68) 



For the reflection coefficients one has to insert from (2.23) according to the case from 
Table (1) with — •>/ {e^^{—u!) — l)u!^ + q^ and the permittivity 

e^^-u) = 1 - . (4.69) 

uiyui — Z7) 

in place of (2.30) accounting for (3.14). 

Now we make in eq. (4.68) the substitution q — > uq. We get a factor u;^ in front. 
The remaining integration remains finite in the limit a; — )■ 0. Thus, the contribution 
from region (a) starts from and is beyond we are interested in. 

In region (b) we have a real 7] and from eq. (2.25) and (4.69), 



x=.L2^^+ 2^ (4.70) 
which is complex for all r]. Thus, we have in (3.31), 

/•oo 

(^(6)(ia;) = / dry 77 In (l - rir2e"'') , (4-71) 

^0 

where we used the convention introduced with eq. (4.9). Here, in distinction from the 
previous cases, the whole integration region contributes. It turns out that it is not 
possible to obtain an expansion, for a; ^ 0, by simple substitution and expansion of 
the integrand. 

First we consider case 1 as defined in Table (1). Here the function (4.71) takes the 
specific form 

/■oo 

V{i,b){ioo) = / ^7777 In (1 + rxEe-'') (4.72) 

with 

rxE = \ (4.73) 

For the expansion we use Appendix A. The integral in (4.72) is the same as in (A.l) 
with the substitutions 

ro^-1, n-'^ul^^. (4.74) 

Its expansion for small Q is given by eq. (A. 3) with the coefficients (A. 19). From the 
leading order term of the expansion (A. 3) we get for (4.72) in leading order in u 
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After inserting into (3.11) we get finally 

= ^ (^E - 1 + ^ In u + 0{u'/'). (4.76) 

Since there is no contribution from region (a) to this order, the above is the complete 
result in leading order. This is just the case involving a logarithmic contribution shown 
in the expansions (3.18) and (3.21). Comparing these with (4.76), we can identify the 
coefficients, 

- % (4-77) 

which, after insertion into (3.21), give the leading order terms in the low temperature 
expansion proportional to and In T. 

Next we consider case 2 from Table 1. Here we have 

poo 

<P{2,b}{ii^) = dr]r]ln{l- TteC"'') , (4.78) 

^0 

with tte given by eq. (4.73). In this case, for the leading order, it is possible to 



make the substitution q — >■ q in the integral. After that the integrand can be 

expanded, 

^ + ^1 Jo (r? + V V + ^) 

The remaining integration is simple and after inserting into (3.11) the result is 

$(2)(u;) = '^^^^^^uj + 0{co'/% (4.80) 

remembering that there is no contribution from region (a). 

The corresponding contribution to the free energy (the first term in eq. (4.83) 
below) was obtained in [19]. In [20], see also [21, 22], also the next order, which is 
proportional to cj"'/^, was considered. The coefficient in front was calculated in some 
approximation. In fact, this coefficient can be obtained exactly using the same method 
as above for the case 1. For this we rewrite eq. (4.78) in the form 

/•oo rco 

n^,b) iico) = dr] T] In {1- rTEe"''/') + / dr] r] In {l + rxEe-"/') , (4.81) 
JO Jo 

reducing the problem, in this way, to the previous case. In both integrals we make 
the substitution r] ^ 2r]. We get a factor of 4 in front of the integral and have to 



34 



substitute ft— > Q,/2 when applying eq. (A.l). Now we use the expansion (A. 3) from 

the Appendix A for the two integrals with Tq = +1 for the first and Tq = — 1 for 
the second, together with the substitution of Q according to eq. (4.74). We restrict 
ourselves to the next-to-leading order (for higher orders one would need to calculate 
the corresponding contributions from region (a) too) and obtain from a2 and in 
(A.19) and (A.21) 

= ^^f- + ^^-"^ + 0(-^). (4.82) 

The leading order coincides, of course, with the result from the integration in (4.79). 
Using eq. (3.19), the corresponding contribution, which results from the TE polariza- 
tion, to the free energy reads 

, _TF 21n2-l<^^ , Cr(5/2) oo^ , , 

AtJ^™= P y2 ^ SR^ / ^ p ^ y5/2 ^o(T^), 4.83 

where we restored the dependence on a. The numerical coefficient Cr(5/2) / (9\/27r^/^) ~ 
0.019 is different from the one in [20]. We mention that, using the formulas from 
Appendix A, higher orders can be written down too. These formulas can be used also 
to obtain the higher orders in case 1, 

= ^ 1 + ^ ^) - + f + (4.84) 

The remaining two cases from Table (1), which belong to the TM polarization, are 
easier if restricting to the leading order. Here we have to insert ttm, eq. (2.24), with 
the permittivity (4.69) into (3.31). The direct expansion of the integrand for a; — > is 
possible. Using, for a; — > 0, 

In (1 - rxMe-'') = (^^7^1)^ ^ + " " " ' (^-S^) 



the integration over rj can be carried out and one comes to 

27r27 
p 

The last, case 4, can be treated in complete analogy and we get 



^(3){u) = -^u + ... . (4.86) 



^(4)M = -^^ + ---, (4.87) 

which completes the calcTiIation of the low temperature expansion (3.21) for the Drude 
model with fixed dissipation parameter. 
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4.4.2 Vanishing dissipation pcirameter 

As already mentioned in the Introduction, there is a problem with the Drude model 
for vanishing dissipation parameter 7. The problem has two aspects. First, consider 
the entropy S, related by the usual relation (1.7) with the free energy (2.1), for T — )> 0. 
In case, the relaxation parameter depends on temperature, 'y{T), and it decreases, for 
T ^ 0, not slower than the first power of T, 

-f(T) ~ r° («>!), (4.88) 

T-^O 

there is a linear term in the expansion, 

= + J-iT + . . . , (4.89) 

of the free energy J-'^'^ of the Drude model which causes a non vanishing contribution 
to the entropy for T ^ 0, 

lim S = -Ti. (4.90) 

T-5>0 

Such non zero J^i constitutes a violation of the third law of thermodynamics (Nernsf 
heat theorem). 

The other aspect appears if considering the free energy of the Drude model for 
vanishing relaxation parameter, 7 — >■ 0, at fixed temperature T. While the permittivity 
(2.27) turned into that of the plasma model, 

e^' ~ (4.91) 

7-5>0 

and so do the reflection coeflicients (2.23), the free energy does not turn into the free 
energy J-'p' of the plasma model, 

jrDr ^ jTPi + jr^j^^ (4 92) 

keeping, in the limit 7 — > 0, an additional contribution, where J-i is just the same as 
in eq. (4.90). This contribution is, from a physical point of view, unsatisfactory since 
a small, or even vanishing, dissipation should have a correspondingly small effect on 
measurable quantities like the force {J^i depends on the separation a). 

It is interesting to mention that the experimental veriflcation of the Casimir force, 
in a number of experiments, favors the plasma model and rules out the Drude model 
[2]. In the experiments the dissipation parameter is small, a typical value is 7= 0.035 
eV while = 9 eV for gold. Now, if assuming that the correct formula describing 
the Casimir force in case of dissipation, is perturbative in 7, it would be clear that the 
dissipation gives only a small addendum wich does not show up in the experiments to 
date. 

There are two ways to derive j^i. The first way uses the representation (2.1) in 
terms of Matsubara frequencies. It is a property of the reflection coefficients (2.24), 

Dr _ Dr _ ^"^(^6)^ " ,^ 93^ 
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with, from eq. (2.31), 



and 



to have a sufficiently different behavior for / 7^ and I — 0. For I — we have — rj 
and £^■'(0) = 00. Hence, for the reflection coefficients 

^?E = 0, r?^ = l (1^0) (4.96) 

holds, while in the plasma model, where we have still £^^(0) = 00, but from (2.28), 



xi = ^u:l + r,\ (4.97) 

a xo 7^ 77, the relations 

r^E-^^^^^p, 4m-1 (l-O) (4.98) 

hold. In this way, the contribution from / = to the Matsubara sum is different for 
both models. It must be mentioned, that in the Matsubara representation, the (/ = 0)- 
contribution to the TE polarization is the only place where this difference shows up. All 
other contributions to the free energy in the Drude model turn into their counterparts 
in the plasma model. Thus, for 7 — > 0, the difference in the free energies between both 
models, 

lim {J^^' - J^^) = J=x T, (4.99) 

results just from the (Z = 0) -contribution in the TE polarization to the free energy of 
the plasma model, see [3]. The contribution from the Drude model vanishes along with 
the reflection coefficient r^^, eq. (4.96). This difference can be calculated easily from 
eq. (2.1), 

1 



with, from (2.2), 



^(0) 



T^^-—^^{^) (4.100) 

/•oo 

/ dkk\^{\-r^{u^,kfe-^). (4.101) 

Here we used, for / = 0, (4.95) and 1] = k which follows from (2.17) and (2.24). Eqs. 
(4.100) and (4.101) were first obtained in [3]. For the reflection coefficient in (4.101) 
wc introduced with 
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Figure 3: The function J^i(Jl) normalized to its value J^i(oo) = — C(3)/47r^, eq. (4.105), 
at infinite Q. 

a separate notation, which will be used several more times below. It results from the 
TE polarization, eq. (4.49), 

n{ujp,rj) = 4^. (4.103) 
The function J^i depends on one variable ft only. We make this fact explicit by writing 



1 f'^ 

(Q) = dkk\u{l-rAQ,kfe~^). (4.104) 

47r^ Jo 



This function can be easily calculated numerically since the integral is rapidly converg- 
ing. A plot is shown in Fig. 3. The asymptotic of this function for large argument, 
which by means of (2.29) corresponds to small skin depths, can be easily obtained from 
expanding the integrand in (4.104) with subsequent integration, 

•^1 = ^ (Cr(3) - 8CR(3)i + 48CR(3)i^ + 16 (-16Cr(3) + Cr(5)) ^ + ■ ■ ■ ■) 

(4.105) 

The asymptotic expansion for small argument is more elaborate and we represent it in 
the Appendix A, eq. (A. 3), where one need to put tq — 1 with coefficients given by 
eq. (A.21). 

It must be mentioned that the above derivation of J^i does not depend on which 
of the two aspects mentioned above, is considered. Moreover, it holds also in the high 
temperature limit which is, in the leading order, for all models given by the zeroth 
Matsubara frequency, 

which vanishes for the TE contribution to the Drude model whereas it is the same in 
both models for the TM polarization. In this way, the TE contribution is missing in 
the Drude model which, for example at large separation, where both polarization give 
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the equal ideal conductor contributions to the plasma model, amounts for as much as 
the half of the free energy. 

The second way to derive T\ starts after the application of the Abel-Plana formula, 
i.e., from eqs. (3.8) and (3.10) with the permittivity £°''(— cu), eq. (4.69), inserted. The 
reflection coefficients are 



^?e(-^,7) = 



TM 



.,2 



tji{u)-i'i) I I y PLO-i-f 



(-^,7) = \ (4.107) 



where we also indicated the dependence on 7 explicitly. We remind, that for 7 — )■ 
these coefficients turn into that of the plasma model. Thus, a difference between both 
models may follow only if this limit and the integrations do not commute. Further 
we mention that both expressions, and J^P', do exist. Both can be obtained and 
coincide if carrying out formally the limit in the integrand. However, as known, this 
is a necessary, but not a sufficient condition. In fact, one needs a uniform convergence 
of the integrand, which is not there. Namely, this limit is not uniform for the TE- 
contribution, 

limr?E(w,7) = n{n), 

7— >-0 

limr?^(cu,7) = 0. (4.108) 

In fact, this is the same behavior as observed in (4.96) in the Matsubara representation. 

Within the representation (3.10) of the free energy, the non-commutativity of the 
limits (4.108) shows up in the TE polarization in region (b), defined in (3.24). In region 
(a), from the integration volume, there is an additional factor a;^ which restores the 
uniformity. Thus we are faced with 

Ar7"(6)TE = ^ ^co/T _ - c.c), (4.109) 

where 

poo 

^M = l rfr^r; In (l-(r?^(-a;, 7))%-"), (4.110) 
and we used (4.103). 

Now, in order to account for (4.108), we split the integration region in (4.109) into 
two parts, 

At^(6)te = ^ + 5, (4.111) 

with 

If' 1 



^=4^ /^ duj^jj^t{^{tuj)-c.c.) (4.112) 
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and 

B = ^,j^ du:-^JJ^^^{^{^u)-c.c.). (4.113) 



In part A we substitute u) cuj, 



7 



1/7 I 



^^4^./o ^^i;3^75^^(^M-^-^-) (4-114) 



with 



= ^ driT] In 1^1 -n {^P\l^^-^'n^ e-^j (4.115) 
which, in fact, does not depend on 7. We used with (4.107) 

r^^{u, u-f) = ri (^^p^^j^, (4.116) 

with ri(r2,r7) defined in (4.102). After this substitution we can put 7 = in the 
integrand and in the upper integration hmit, or tend T ^ assuming a temperature 
dependent 7(T) with hmr^.o = 0, in (4.114) and obtain 

A{^^0) = ^J^ -^i{i;{^u)-c.c.). (4.117) 

In part B we can put 7 = directly in the integrand and in the lower integration limit 
and obtain the corresponding contribution to the plasma model. 

In this way, we re-obtain eq. (4.99) with the right hand side given by A, eq. (4.117). 
The last step is to carry out the integration over u in eq. (4.117). This is possible 
using eq. (B.5) accounting for 99(0) = 0, which holds for the function (4.115), and we 
get 

A{^^0)^J^iT (4.118) 

with J^i given by eq. (4.104). 

In fact, the division (4.111) of At'J^(6)te was observed in [23] in an attempt to pass 
from the Matsubara summation to an integration over real frequencies, just the same 
way as one does applying the Abel-Plana formula. In [23] it was shown what happens 
to the integration path under deformation, ^ = cje~*'^, with /3 changing from /3 = 
to /3 = 7i/2. It was found (see Figs. 5 and 6 in [23]) that this path with increasing 
/3 develops a self- intersection and that, for ^ = 7r/2 in the hmit of 7 — >■ 0, it decays 
into two separate pathes which just correspond to the two contributions, A and B, in 
(4.111). In [23] this contribution was found analytically. 
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4.5 Insulator described by oscillator model 

For an insulator, we use the permittivity given by eq. (2.32) and for the low tempera- 
ture expansion we need to insert 

£i-ui.(_^) + y — ^_ — (4.119) 

into ip{ibj\ eq. (3.31). Like in Section 4.2, with the permittivity (4.119), the free 
energy of an ideal conductor does not follow in any limit. Therefor we use in this 
section also the four cases, introduced with Table 1. 

We start with the special case of no dissipation, 7^ = in (4.119). This is a real 
permittivity. We expand it for small a;. 



£-«"i-(-c^) =£0 + ^ + ... (4.120) 



with 

£o = 1 + E4> -2 = y.% (4-121) 

The absence of an imaginary part in (4.120) allows to use the results obtained in 
Section 4.2 for fixed permittivity and simply to insert (4.120) into (4.31), (4.33) and 
(4.32) with subsequent re-expansion in powers of io. Since the expansion in (4.120) is 
in even powers of a;, the first correction due to (4.120) occurs in the order u;^, and the 
next in order cj^. Denoting the coefficients of the expansion according to eq. (3.18), 
one has to substitute, besides £ ^ £0, in the coefficients 

1 

$(fe),4 ^ $(fe),4 + -2$(fc),2, 

^(fe),5 ^ ^(fe),5 + A*(fe),2, (4.122) 

where (A;) denotes the case according to Table 1. 

Now we allow for dissipation, 7^- 7^ 0. In this case the expansion of the permittivity 

is 

£--i-(-a;) =£o-i7o'-' + --- (4.123) 

with 

70 = (4-124) 

Cross terms, involving corrections from both, (4.120) and (4.123), do not occur in the 
orders which we are interested in. 

We divide the integration in Lp{iuj) into the regions (3.24). In region (a) we use 
the same formula, eq. (4.12) as in Section 4.2. We can go the same way as there 
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and obtain an expansion which starts from uu^ which will be non leading. The leading 
contribution comes from region (b). Here we have to use eq. (3.31), 

/•oo 

(/3(b) (icu) = / dr] r] In (l- rir2e"'') , (4.125) 
Jo 

with the reflection coefficients (2.24) to be inserted according to the cases in Table 1. As 
it would be too elaborate and of less use to get higher orders of the expansion for small 
u, we restrict ourselves here to the leading order, which is u, and to the logarithmic 
contribution in order co^, i.e., the contribution being proportional to InuJUJ^, which 
is easy to obtain. We do not calculate the contributions proportional to uu^ without 
logarithm. The simple reason for this restriction is in the structure of the integrand in 
eq. (4.125), where an expansion in powers of u delivers powers of r] in the denominator 
which make the integration over rj divergent. Also, a substitution rj ^ rju is not helpful 
since it produces a problem on the upper interaction boundary. The leading order can 
be obtained from expanding the logarithm in (4.125) in powers of ou, 

In (1 - nrse-'') = f,{ri)u + h{r])u:^ + . . . . (4.126) 

We consider here odd powers of uj only since the even powers will drop out later 
from the imaginary part. The function f\{r]) is integrable, whereas the function f3{r]) 
behaves as 

fM = H + 0{l) (4.127) 

for — )■ 0. This behavior results in a logarithmic singularity when inserted into (4.125). 
As a consequence, the expansion of the function <^{uj) has a logarithmic contribution, 

$(fe)(a;) = + ($(fc),3 + ^[k),^ uj^ + 0{uj^) (4.128) 

where we used the notation of eq. (3.18). The coefficients $(fc),i can be obtained by 
simple integration from /i(r/) for a specific case according to Table 1, 

/•oo 

$(fe),i= / driTjiiMrj)- c.c). (4.129) 
Jo 

For the coefficient ^{k),3 we split the integration, 

/•oo ru! poo 

/ drjrjiifsirj) -C.C.) ^ drj rj i{f 3(7}) - c.c.) + drj rj i{f 3(7]) - c.c.) . (4.130) 

Jo Jo Jlo 

The leading contribution comes from the lower integration limit in the second term of 
the right side and we get 

*(fc),3 = -^(/30-c.c.). (4.131) 

As said already, we do not calculate $(fc),3. The results are shown in Table 3. These 
contributions are all due to dissipation. In order to get the corresponding results for 
the electromagnetic case, the polarizations must be added, cases 1 and 3 for an ideal 
conductor in front of an insulator and cases 2 and 4 for two insulators. 
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case 



^(fe),3 



2 



1 












4 



3 



7o 



7o 




Table 3: The coefficients $(fe),i and ^{k),3 calculated in Section 4.5 for an insulator with 
dissipation. 

4.6 The case of dc conductivity 

In this subsection we investigate the contribution of a dc conductivity to the free energy. 
For imaginary frequency, the permittivity is given by 



assuming Eq ^ 1. The behavior of this model is quite similar to the Drudc model 
considered in Section 4.4. Here, it is the conductivity a, which describes the dissipation 
of energy. The counterpart, in which the free energy corresponding to (4.132) does not 
turn, is the insulator considered in the preceding subsection in place of the plasma 
model in Section 4.4. In the following subsections we consider separately the cases of 
non vanishing and of vanishing conductivity a. 

4.6.1 Fixed conductivity a 

The case of a fixed conductivity a, not depending on temperature, is in fact nonphysical 

and of rather academic interest. We include it here since the calculations go in parallel 
to a fixed dissipation parameter 7 in the Drude model, only with the role of the two 
polarizations interchanged. In both cases the permittivity behaves the same way for 
small frequency, it has ou in the denominator to the first power. The difference between 
both is in the behavior for large frequency, where e'^^{uj) — >■ £0 with £0 7^ 1, whereas 
e^^{uj) — >■ 1 in this limit. 

For the calculation we use the division into cases, given by Table 1, and into regions 
(3.24). We are interested in the leading order only which is co, including culncj. The 
region (a) gives contributions starting from cu^, which we do not consider. Following 
(3.31), in region (b) we have to consider 



with ri and r2 to be substituted according to Table 1 by one of the refiection coefficients 



Ana 



(4.132) 




(4.133) 
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from (2.24) with the permittivity 

e'"^{-cu) = eo + ^. (4.134) 

lU! 

and X from eq. (2.25). 

The further considerations are in parallel to the Drude model, Section 4.4., and 
we do not repeat all discussions. For the case 1, for the function (p(^i-^{iu!), we use eq. 
(A.l), where we have to substitute 

ro^-1, n'^ ^ Aniau - {eo-l)uj'^. (4.135) 

From the leading order coefficient from eq. (A. 19) and taking the imaginary part 
according to eq. (3.11), we get 

$(1) (uj) = 27r(T (^-1 + 7e + ^ ln(167r(7a;)^ cu + 0{uj^''^) , (4.136) 

which has a logarithmic contribution. 

Case 2 is easier. After substituting rj \/ojr] in (4.133), one can expand the 
integrand, 

rd„l.-ip^^^ + 0(.^V.) (4.137) 

with subsequent integration. This results in 

$(2)(c^) = 27r(T(ln2 - \)uj + 0{uj^''^). (4.138) 

In the leading order, considered here, the results are in parallel to the Drude case, 
considered in Section 4.4.1. This can be seen comparing the substitutions (4.74) and 

(4.135). To leading order these are connected by the substitution — > Aira. This holds 

also for the next-to-leading order, which is proportional to a;^/^ and which results from 
the coefficient 03 from the Appendix. In this way, the formulas (4.82) and (4.84) 
from the Drude model translate by the above transformation into the corresponding 
formulas for fixed dc conductivity. 

In the cases 3 and 4, i.e., for the TM polarization, one expands the integrand 
directly for small ui, 

m>('^) = fd, ^^.^^l_^^ +0(.\ (4.139) 



which results in 
In case 4 the result is just twice. 



^{3)H = + 0{uj^). (4.140) 

6(7 



$(4)(cu) = -^cu + 0{uj^). (4.141) 

These results are similar to that of the Dmdc model. Section 4.4.1, and can be obtained 
from (4.86) and (4.87) by the above substitution. 
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4.6.2 Vanishing conductivity a 

In this subsection we consider the free energy for vanishing dissipation parameter, 
(7 — > 0. As before in the Drude model there are two aspects. The first is the entropy 
eq. (1.7), having a residual value 

lim S = -J'f (4.142) 

T-5>0 

similar to (4.90) for vanishing temperature in case a goes to zero not slower than the 
first power of the temperature such that 

lim ^ = (4.143) 

holds. The second aspect is that the free energy with dc conductivity, J-''^^, for o" — )■ 
at fixed temperature, does not reproduce that of the insulator J^^^*^'- considered in 
Section 4.5., 

jrdc ^ -pinsul. ^ -pa J. ^ (4.144) 

where J-'f is the same as in (4.142). The comments made on this in the Drude model 
apply here too. 

There are two ways to show the above statement and to calculate J^f . The first way 
uses the Matsubara representation. The reflection coefficients are given by eq. (2.24) 
with (2.25) and the permittivity (4.132), 

V + V £o^f + 477(7^, + 772 

(so + ^) - vW+4^^^6T^ 
4^m(6,<^) - i , (4.145) 



and the corresponding coefficients for the insulator are given by the above for a — 0, 



insul. 
TM 



(a) = -'"-fgp! (4.146) 



for all /. Now, there is a difference between both for zero Matsubara frequency, / = 0. 
For dc conductivity we have 

4^e(0,(t)-0, 4l,{Q,a) = l, (4.147) 

whereas, for the insulator, 

rVr-(0) = 0, 4°r(0) = J^, (4.148) 
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hold. As compared to the Drude model, in this case the TM polarization delivers the 
difference. In all other contributions the limit cr — > is smooth and we get 



(7-)-0 

with 



with 



lim (J^" - ) = J=IT (4.149) 

■f = — / d^^ (In (1 - e-") - In (1 - r^e"")) (4.150) 
^'^ Jo 



.0^1^. (4.151) 



The integration can be carried out delivering 

^1-^2 (Cr(3) - Li3(ro')) (4.152) 

(see [24]), which comes in place of (4.104) in the Drude model. 

The second way to derive J-'f starts from representation (3.8) with (3.11) and the 
permittivity (4.132). The reflection coefficients are 



rTE(-a;,(j) = 



,2 



4m = -7^ ; =• 4.153 

(^0 + ^)V+ V^ooj' - ^TTiau + T}^ 

The coefficient of the TM polarization does not have a limit for cr — t- which would 
be uniform in u. Again, the interesting contribution comes from region (b). 



which we split also, 
with 

and 



— j^ d77 77 ln(l-(4U-a;,<7))%-''), (4.154) 

At^tm = ^ + ^, (4.155) 
I 1 



The part B turns, for u — > into the corresponding one of the insulator. In A we 
substitute uj uja, 

A = ^ f^''^ ^^^j^ z(^(z^) - C.C.) (4.158) 
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with 

drir) In \l-ri{eo + —,v] e"" 1 . (4.159) 



roo 



The function ri(Q, rf) is given by eq. (4.102) as before and '^{iu) is, in fact, independent 
on cr. 

Now we can tend cr ^ in (4.158) and get 

A(a^O) = — — tiMtuj) - c.c.) (4.160) 
47r^ Jo ^ 

The last step is to apply eq. (B.5) with (p{iu) which results in 

A{(T ^ 0) = J^^T (4.161) 
and finishes the second way to derive J^f . 

4.7 Hydrodynamic model for graphene 

The hydrodynamic model is described by the reflection coefficients (2.35). Like in the 
plasma model in Section 4.3, we allow for different plasma frequencies, LOpi {i — 1,2), 
in the two interfaces. The corresponding reflection coefficients are 

-1 1 , , 

rxEi = ^_ r^Mi = ^ 2 ■ (4.162) 

For the calculation of the low temperature expansion we use again the division into 
the two regions (3.24). 

4.7.1 Region (a) 

Here we have a real q and the function (3.30) is 



/•ui 

/ dq q hi {1- rirae"^*) , (4.163) 
Jo 



where for one needs to insert r^Ei or ttmi according to the polarization. In this 
model, the calculation of the expansion of for small u is not such easy as for the 
plasma model in Section 4.3.1, for instance, the reflection coefficients are not pure 
phase factors. 

For the TE polarization, however, the expansion can be obtained easily since an 
expansion of the integrand in powers of cu does not cause problems with the conver- 
gence of the subsequent integration over q. In this way one can obtain the expansion 
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machinized and insert into (3.11), and carrying out the integration the same way one 
obtains 

*(a)TE('^j = "o"^ + I o 7 , , \ 1 ^ 

Z \6 OUJpiUJp2\UJpiU!p2 + UJpi + UJp2) J 

+0(u;^), (4.164) 

where we restricted ourselves to display the first two contributions only. 

For the TM polarization such a simple expansion does not work since an expansion 
in powers of cu would produce inverse powers of q which make the integration divergent. 
We proceed as follows. We split the functions </^(a)TM(^'^) into two parts, 

<^(a)TM(iCc') + (4.165) 

with 

A= dq q In {1- rTMirTM2) (4.166) 

^0 



and 

B 



f dqqlnfl-- ^ ). (4.167) 

Jo \ ('^TMl?^TM2) - 1/ 

The integration in A can be carried out explicitly with subsequent expansion in powers 
of a;. In B, the integrand can be expanded up to o;^ with convergent subsequent inte- 
gration. Adding the results from A and B and using (3.11) one obtains the expansion 



*(a)TM('^) = -^^'^ + 



^3 UJpi Up2 {u!plUp2Y/'^y/Upi + Up2 ^ 



TT (^^1 + 2u:l^U^2 + i^lii^l2 + 2c^pl^ |2 + ^p2) ^.4 ^ q(^5)^ 



'^^ll^Ui^pl + ^p2f 

which completes the contributions from region (a). 



4.7.2 Region (b) 

In this reagin we have a real 7] and, inserting q = it] in (4.162), real reflection coef- 
ficients, which must be inserted into (4.51). Thus, just as in the plasma model, we 
have a seemingly real function 93(6) (io;), eq. (3.31). It can have an imaginary part 
only from the logarithm when its argument changes sign. Indeed, the hydrodynamic 
model, in parallel to the plasma model, is known to have surface plasmons in the TM 
polarization. In order to account for them we rewrite eq. (3.31), 

poo 

(P{b){iu!) = / drjrj (-ln(rTMi)~^ - ln(rTM2)~^ + In ((rTMi)~^(rTM2)~^ - e"'')) , 

(4.169) 
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removing the poles from the arguments of the logarithms and separating the contri- 
butions from the plasmons on the interfaces taken alone, 77singie,i — 1)2), which are 
solution of 

(rxMi)-' = (4.170) 
and the two plasmons r^symm. and r/antisy. which are solutions of the equation 

(rTMi)"'(rTM2)"' = e-''. (4.171) 

We took the notations in accordance with the symmetry the corresponding wave func- 
tions have in case of equal plasma frequencies. 

With (4.162), the solutions of equation (4.170) are obviously 

Vsmg\e,i = ■ (4.172) 

The equations (4.171) are transcendental ones like in the case of the plasma model. 
Section 4.3.2. There is only a small difference to the case of the plasma model as here 
the solution r^antisy. exists also for arbitrarily small cu. 

The solution ?7symm. can be obtained from the equation (4.171). We solve this 
equation for rj occurring in the left side, which amounts in solving a quadratic equation. 
With the appropriate one of the two solutions we rewrite the equation in the form 



l-yi-2//(l-e-'') 



2 



rj = ^ ,^ , , (4.173) 

//(1-e-^) ujpi+ujp2 

where we introduced the notation 

(4.174) 

{ujpi + ujp2y 

This equation can be soled by iteration, starting from inserting 77 = in the right side. 
After re-expanding in powers of ou, the solution is 

= + I ( -^)\ ( -^)' + 0(.'). (4.175) 

c^pi+i^p2 2\u;pi + u;p2j 4 \ujpi + u}p2 J 

In order to obtain the solution r/antisy.; one needs to take the other solution of the 
quadratic equation mentioned above and one needs to rewrite the equation in the form 



= |^^l + Vl-2^^{l-e-v) u;p,+up2^_ 
Iteration, starting from inserting 77 = in the right hand side, gives 



UJpl+UJp2 , h^pl+(^p2 
^antisy. - \ -\ a U 

*' <^pi<^p2 4 ly u;piu;p2 
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Using these solutions, we return to eq. (4.169). The analytic continuation can be 
carried out starting from large rj, where the logarithms are real. The arguments of 

the logarithms change sign whenever, for decreasing 77, these pass a zero, and the 
logarithms acquire an addendum ivr. The argument of the last logarithm in (4.169) 
has two zeros, r^antisy. and ?7symm., whereby the inequality 

?7symm. < ?7antisy. (4.178) 

holds, which can be verified for small u with the leading terms in (4.175) and (4.177). 
Since there is no pole in between, the logarithm acquires an addendum of —iir in 
passing rjsymm., i-e., the second zero. In this way we get from (4.169), and accounting 
for (3.11), 

(/•^/single,! /■»?single,2 A'»)antisy. \ 

- / drirj- / driri+ / dr]r]\. (4.179) 

Jo Jo Jrisyrara. J 

Carrying out the integrations and re-expanding for small a;, we get 

TT / 1 1 \ 2 , T^i^ll+^h) 3 



2 \Upi U!p2 



$(,)(a;) = _ + _ + , ^ (4.180) 



TT (91a;^i + 184cu3ia;p2 + 294a;2^a;2^ ^ \Mupx'42 + 91^p2) 4 , ^/ 5x 

96(^;i(^;2('^pl +'^p2)^ 

This is the complete contribution from region (b) since the TE polarization does not 
have surface plasmons. 

4.7.3 The low frequency expansion 

Here we collect the results from the preceding two subsections. For the TE polarization 
we have only the contribution from region (a), which is given by eq. (4.164) 

^TE(t^) = $(a)TE('^)- (4.181) 

Special cases are equal plasma frequencies, 

$te((^), = _7r 2 2(a;^ + 3) 3 

^ ^kpi=-p2=-P 2 3a;p(wp + 2) ^ ' 
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30a;3(a;p + 2)3 
and one interface ideal conducting, 

$teMi = --u^ + ^^^f^^u^ (4.183) 

< + 8a;,^ + 2Aul + 24a;, + 8 
60.|KH-1)3 " 

For the TM polarization we have to add (4.168) and (4.180), 



2 V i^pi i^p2 



3 ^Pi 4V2^a;3,a;32(a;,i+a;,2) 
(a;2i+a;2_^)arctan/' /J££l^^ \ 



LUp 1 +LJp 2 

(<^pl<^p2)^/^V'^pl + 0Jp2 



oj'^ (4.184) 



TT 



.. (43a;,^ + 88a;3,a;j,2 + ^A&ojI^uI^ + 88a;pia;32 + 430;^%) 

^ ^ n« 2 2 ^ ; ^ ^ +0{UJ). 

Special cases are, again, equal plasma frequencies. 



Up 



TT ^ , 2 \ 2 /I 2 TT \/2arctan^^^^ _^ 



$tm(u;), _ = -- 1 + — + - + + ry^^ ^ a; 

^\.p,-^p,-^p 2 V '^p/ \^3 Up 4^3/2 ^3/2 I 

1277r / -16a;3 + 47a;^ + 402a;, + 480 
^^'^ + I 96a;3(a;, + 2)^ ^^"^^^^ 



5^ 5arctan(^^y . 70^^ 



512a;p'/' 32x/2a;,^/' / 11520a;| 
and one interface ideal conducting. 



TT 



^ , , ,. 2 / - - arctan(y^;) , 

/|.pi=oo,.p2=.p 2 l3 4V2a;p/^ ul'^ 



437r 4 / 57r -8a;3 + SQw^ ^ 113cUp + 64 
+ 96^"" + 1^"2567^ 96a;3(a;, + If 
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5 arctan (^/u^) 



512a;. 



5/2 



IQtt 
5760a| 



cu6 + 0(cu^). 



(4.186) 



Finally we add the polarizations and get the expansion for the electromagnetic case, 



1 1 

+ 



p2 



TTU! 



(4.187) 



+ 



Su!piUp2{u!plU!p2 + <^pl + ^p2) 

^(,,2 ,^,2\ (ujI,+ ujIo) arctan ( f^I^ 

T^y^pi+^^p2) ^ P'^ P^-' yy ojpi+ujp'i 



+ 



+ 



A^^ul,ul2{^p,+Up2) (^Pl^p2)3/V^pl+^P2 

TT (43a;^i + SSojl^Up2 + 246cc;2ia;22 + 88a;pia;3^ + AZu^^) 



I 

a;^ + 0(a;^). 



Special cases are considered in the next section. The coefficients of these expansions, 
identified according to eq. (3.18), must be inserted into eq. (3.21) to obtain the low 
temperature expansion of the free energy in the corresponding case. 



5 Compilations of results 

In this section we represent the results for the low temperature expansions of the 
free energy for the specific models defined in Section 2. The relation between the low 
frequency expansion, eq. (3.18), of a function $(0;), introduced in eqs. (1-14) or (3.10), 
and calculated in Section 4 for the specific models, and the low temperature expansion 
of the free energy is given by eq. (3.21). Also, here we restored the dependence on 
the widths a of the gap according to the rules given in Section 2, eqs. (2.5)-(2.7). 
However, we did not restore the dependence on the constants c and k-Q since that 
would overload the formulas and refer to the corresponding rules in Section 2, eqs. (2.8) 
and (2.9). Also, in this section we do not consider the polarizations separately. All 
formulas are written for the electromagnetic case, i.e., after adding the contributions 
from both polarizations. We follow the order of the models used before. 



5.1 Ideal conductor 

We start with the case of ideally conducting interfaces considered in Subsection 4.1. 
Prom eq. (4.8) we get 

^^id.cond. ^ _ Cr(3) y3 ^ ^ y4 ^ _ _ 

27r 45 

where the dots denote exponentially decreasing contributions. As well known, the 
leading order is independent from a and it does not contribute to the force. 
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5.2 Fixed permittivity 



Here we consider dielectric medium with permittivity. For dielectrics behind both 
interfaces we get from eq. (4.32) 

£ + 1 An ^ ^ ^ ^90 

47r 

+{e -l){e- 1)(24£=^ - 36s' + 29s - 6)-— 
+ (£ _ 1)3(^2 ^ a'r + 0{ T') . (5.2) 

In case of one interface ideally conducting we get from eq. (4.33) 

£ + 1 Stt ^ ^ ^ M5 

_(^_l)2£23CR(5)^2y5 
TT 

+ (£ - l)2£3/2(24£3 _ 36£2 + 5£ + 10) ^ a3 

+(£ _ 1)3^2 SOCr^TK t^ + o{T'). (5.3) 

Like in the case of ideal conductors, the leading order does not contribute to the force. 
The first contribution to the force comes from the fourth order in T. In eq. (5.2), the 

coefficient in front of is, up to a factor, C4 in eq. (12.97) in [1] and the corresponding 
coefficient in eq. (5.3) is, up to a factor, K4 in eq. (15.17) in [1]. The contributions 
of orders and in eq. (5.2) were obtained in [24, 25]. In [24] the coefficient C4 is 
also calculated for dissimilar plates. 



5.3 The plasma model 

We start with the case of media with different plasma frequencies, Upi and u)p2, behind 
the interfaces. Using eq. (4.65) we get the expansion 



a \0Opi ujp2 J / 27r \ a \0Opi u}p2 J J 45 



2 \ujpi ' Up2j a oj^2 
'up2 _}_\ \ 2Cr(5) ^5 

+ {^ + ^^^T^ + 0{r). (5.4) 



^12) 105 
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It is represented in a form that the behavior for large plasma frequencies can be seen 
easily. For both tending to infinity the ideal conductor case is recovered. The formula 
simplifies if the plasma frequencies behind the interfaces are equal, oupi — ujp2 — oup, 



acjpj 27r \ acdpj 45 



+ fi---4i)^r= + -|!Lr. + o(n (5.5) 

\ aoOp a^Up J UpTi 1050;^ 

The case of one interface ideally conducting can be obtained from tending one of the 
plasma frequencies to infinity, 



2n \ ttujpj 45 

+ {l- — -4-y-^T^ + ^^T^ + 0{T\ (5.6) 

\ aujp a^cOpJ cUpTi 105cj^ 

and the index on the other frequency was dropped. The contributions of order to 
were obtained in [26]. These contributions, with account for (2.29), reproduce eq. 
(14.14) in [1], which was obtained as an expansion in powers of 6 = l/up. In general, 
the free energy for the plasma model has an asymptotic expansion in both, T and S. 



5.4 The Drude model 



Here we consider the case of a fixed dissipation parameter 7 (for vanishing 7 see Section 
4.4.2.) and restrict ourselves to the lowest order which is T^, including T^lnT. For 
medium behind both interfaces, adding the contributions from eqs. (4.80) and (4.87) 
we get 



-^^(21n2-l)- 



7 



27r^7 



^ + o(r^/^). 



(5.7) 



in agreement with [19]. In case, one interface is ideally conducting, we have to add the 
contributions from eqs. (4.77) and (4.86) and we get 



967 



, . s 1 8au„ 

7E + Ck(-l)-l + -ln-^ 



7r^7 



, )2 
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which has the logarithmic contribution. 



ln(47rar) T'^^0{T^>'^), 
(5.8) 



5.5 Insulator described by oscillator model 

For an insulator described by the oscillator model, the expansions follow from Subsec- 
tion 4.5. In case without dissipation, from eq. (4.122), we get 

AtT = -—-——-T^{^,-l)[e,^e, -2j-aT 
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+(£o - 1) ((^0 + 1)(24£3 - 36el + 29eo - 6) 

(5.9) 



and for one interface ideally conducting, 

ArT = -i£^%e)T3+(.f(eo-2) + l)^aT^ 
£o + 1 Stt V ' / 45 

1'° + (.0 + l)2aHV ^ 
+(£o - 1) (£o^'(24£3 - 36el + 5eo + 10) 

-'-^^)^a'T^ + 0(T^). (5.10) 
£o + 1 Sa^cu^ j 2835 ^ V ; V ^ 

These formulas differ from (5.2) and (5.3) only by the substitution e — )■ £o and the 
contributions depending on ujq, i.e., starting from the order T^, which is quite trivial, 
but was not considered before. The first two terms in (5.9) and (5.10) were obtained 
in [24, 25]. 

In case of dissipation we consider the leading order, which is T^, and the order 
T^lnT as another example for logarithmic contributions in the low temperature ex- 
pansion of the free energy. We did not calculate the order (without logarithm) . For 
medium behind both interfaces we get from Table 3 adding the second and the fourth 
lines and using eq. (3.19) 

(5.11) 

and, for one interface ideally conducting, from adding the first and the third lines, 

^0 - 1\ 70 ^2 (^0 - 1)(£0 + 3)7o TT^ , „s „4 , „.„4s 



(6.12) 

where we accounted for that 70, eq. (4.124), has dimension of one inverse frequency. 
The contributions from frequency dependence and dissipation are additive in the con- 
sidered orders of the expansion. The logarithmic term contributes to the force, whereas 
contributions proportional to 70?"^ (without logarithm) do not contribute due to a can- 
celation after restoration of the dependence on a. Cross terms, depending on both, ujq 
and 7o, appear in higher orders, starting from T^, only. In the above expansion, in eq. 
(5.11), the order was obtained in [27]. 
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5.6 The case of dc conductivity 

We present here the results for a fixed conductivity cr, which is rather nonphysical, for 
completeness (for vanishing a see Section 4.6.2.). The low temperature expansion of 
the free energy follows, for medium behind both interfaces, from adding eqs. (4.138) 
and (4.141), 

7r(j(ln2 — 1) TT 



12 



288aV 



(5.13) 



The contribution with a in the denominator results from the TM polarization. It was 
found in [28]. For one interface ideally conducting we get from eqs. (4.136) and (4.140), 



At J" 



ff(7. + CU-l)-l + iM32.a.))-5^ 

+^ln(47rar) T^ + 0{T''^^), 



rp2 



(5.14) 



which has a logarithmic contribution like eq. (5.8) in the Drude model. Both these 
expressions follow from (5.7) and (5.8) by the substitution — >■ Ana. 



5.7 Hydrodynamic model for graphene 

Here we have again two different plasma frequencies in the interfaces. From eq. (4.187) 
we get the expansion 
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which is similar to eq. (5.4) in the plasma model. The special case of equal plasma 
frequencies is 
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Again we see that for cUp ^ oo the ideal conductor case is recovered. Finally we note 
the case if one interface is ideally conducting, 
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It is interesting to mention that the order does not depend on the width a. 



6 Conclusions 

In the forgoing sections we considered the asymptotic expansion of the free energy as 
given by the Lifshitz formula for low temperature. Wc used the Abel-Plana formula and 
exploited, in a unified treatment, the analytic properties of the refiection coefficients. 
Of course, we reproduced the known results, in a number of cases we were able to go 
beyond. 

For instance, in the plasma model, it is possible to get arbitrarily many terms of the 
expansion in a systematic way. In other models, more elaborate methods are necessary 
to calculate higher orders in the asymptotic expansions of some integrals. The basic 
tool, used in this paper, is eq. (1.14) together with eqs. (1.15) and (1.16). Together 
with the division (3.24) into regions, this allows for some more insight into the structure 
of the free energy at low temperature. So, there are from region (a) contributions with, 
besides real frequency real wave numbers q and ^3, which correspond to scattering 
states in the spectrum. From region (b) there are contributions with real frequency and 
imaginary wave number, q = ir\ (inside the gap), which correspond to surface plasmons. 
In this way,the low temperature expansion is determined by the low lying excitations 
in the spectrum. In this sense, the expansions for the specific models, considered in 
this paper, can be viewed as illustration to the general formulas in Chapter 5.2 in [1]. 
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It must be mentioned that this picture becomes significantly more comphcated if 
inchiding dissipation, which brings its own contribution to the imaginary part and 
contributions from all frequencies, including high ones, play a role. An example is the 
Drude model in Section 4.11 with fixed dissipation parameter 7, where in case 1, in 
order to obtain the expansions (4.75), eq. (A.l) from Appendix A was used. To that 
formula, arbitrary high wave numbers contribute. 

In this paper, in the cases with dissipation, i.e., the Drude model and in the case 
with dc conductivity, we re-derived the contributions violating the third law of ther- 
modynamics. These models are at once non perturbative in the respective dissipation 
parameter in the sense that, for vanishing dissipation parameter a contribution is left 
in addition to what one calculates starting from no dissipation. These contributions 
are derived in the original way starting from the Matsubara representation, and in an 
alternative way starting from representation (1.14) in terms of real frequencies. Since 
we do not enter here the discussion on the applicability of these models together with 
the Lifshitz formula, we restrict ourselves to the following remark. As known from 
dielectrics with fixed permittivity, the low temperature expansions and the limit of the 
corresponding parameter, like £ — >■ 00, do not commute. The same happens in the case 
with the dissipation parameters, 7 or o", tending to zero for the corresponding 'danger- 
ous' polarization. So, for instance, the TE cases in the Drude model. Section 4.4.1, cqs. 
(4.76) and (4.80), and the TM cases in the dc model in Section 4.6, eqs. (4.140) and 
(4.141), have the dissipation parameter in the denominator. In the correspondingly 
other polarization this parameter is in the numerator. The latter holds also for 7 in 
Section 4.5, Table 3. 

We included also the hydrodynamic model for graphene. Section 4.7. This model 
is very much like the plasma model, which is quite natural since both can be derived 
from a charged fiuid, one time confined to a half space, the other time to a plane. 
So, in the low temperature expansion, we see similar structures, especially concerning 
the role of the surface plasmons. As for the applicability of these results to graphene 
we mention that its low frequency properties are better described by the Dirac model 
which is not included at this time. 
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A Asymptotic expansion of an integral 

In this Appendix we calculate the asymptotic expansion of an integral appearing in 
Sections 4, eqs. (4.72) and (4.104), and in Section 6, eq. (4.133) for k — 1. It can be 
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written in the form 

poo 

= / ^7777 In (1 - rori(Q)e-'') , (A.l) 
Jo 

where Tq is a number and 

nm - (A.2) 

was introduced in eq. (4.102). 

We are interested in an expansion for Q — > up to some power , 

no /0\^ 

7(Q) = J]a,(ro) - (A.3) 
1=0 ^ ^ 

A direct expansion of the integrand is not possible since it would be in powers of 
resulting in a singular integration at — )■ 0. The same way an expansion after the 
substitution 1] ^ rju is not possible since it would result in singularities for 77 — >■ 00. 
We proceed as follows. First we make a substitution, 

V = ^{t-t-'), (A.4) 

in (A.l), 

I{n) = tLj^ rftt(l-r^)ln(l + ^e-t(*-V*)). (A.5) 
Next we assume |ro| < 1 and expand the logarithm, 

-rn 



/(0) = J]^^(0) (A.6) 
with 

A^{n) = (^^y l^^dtt'-''{l-t-')e-f('-'/'l (A.7) 

In this representation, for small Q, the main contribution comes from large t. Thus we 
can expand the exponential 

\ / fe=0 ^ ^ 

which gives k additional powers of Q in excess of the factor in front of the integral 
and it is sufficient to take ko — riQ — 2. We insert this expansion into (A.7), 

Asm = (§) ' f dt t^-^^ (1 - t-^) e-f^ E ^ (f ) ' + 0(0-+^). (A.9) 
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We split this integral into two by adding and subtracting the integral from t — till 

As{n) = Ps{n) - Qs{n) (a.io) 

with 
and 

Qsm = ) ' I' dt (1 - t-) e--* i: ^ (f ) ' + (A.12) 

Making the analytic continuation to sufficiently small s, the integral in (A. 11) can be 
done, 

- E ^ (f ) (r(2 - 2. - - ) r(-2 - 2. - j . 

(A,13) 

The split in eq. (A.IO) is possible for s < Eq. (A. 13) provides the continuation 
to the whole s-plane. The poles appearing in s = ^, = 0, 1, . . . , will be com- 
pensated by corresponding poles in Qg{fl), leaving logarithmic contributions behind. 
These poles appear for s = 0, 1, . . . , Sq with sq — [^]- All relevant powers of Q come 
in Ps{^) from s < Sq. 

In Qs{fl), (A.12), because of t < 1, we can expand the remaining exponential, 

Qs^ = (f )' dtt^-'^ (1 - t-') E ^ (f - t-')^ + 0(0-+^), (A.14) 

and the integration becomes simple allowing for an easy machinized treatment. With 
these representations for Ps{fl) and Qs{fi), we return to eq. (A. 6). 
We split the sum, 

7(Q) =^ + S + 0(Q"°+^), (A.15) 
into ^ 

A = f2 (Psi^) - (A.16) 

.=1 * 

and 

^=E^Q^(^)- (A-17) 

so+l 

The sum in A is a finite one and that in B is over a rational function of s. For any given 
no, these can be easily calculated machinized. In doing so one can generate expression 
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for the coefficients ai(ro) in (A. 3). The first two coefficients are zero, ao(ro) = and 
O'i{fo) — 0. The next two are 

a2(ro) = ^ h — — ln(l -ro) + roln-, 

«3(?^o) = o ^ 372 arctanh^ro + 7e- (A. 18) 

Higher orders can be easily generated. For special cases for tq the expressions simplify. 
For ro = — 1 we get 



a^{-l) = -log(21])-7E + l, a3(-l) 



8 
3' 



128 
2 ' 



a4(-l) = 21og(f)) + log(2) + 27E-3, a^{-l) = -f, (A.19) 

ae(-l) = -|log(f)-M-41og(2)-«f + g, a,(-l) 
for ro = we get 

a2(0) = log(f)+7E+i, 03(0) = 



32 
'15' 



04(0) = -log(Q)+log(2)-7E + i|, 05(0) = Iff, (A.20) 

512 
' 4725 ' 



06(0) = gig (I201og(f) + 1207E-I9l), 07(0) = 
and for ro = 1 it is 

02(1) = log(§)+7E, 03(1) = - 



04(1) = log(2), 05(1) = -|, (A.21) 

06(1) = -iM + iM^ ^^(1) ^ _^ 



225' 



where 7e is Euler's constant. 



B The (j-integration in eqs. (4.117) and (4.160) 

We consider the integral 

7= r ^i^^(iuj)-(p{-iuj)) (B.l) 
Jo ^ 

as it appears in eqs. (4.117) and (4.160). Since the functions Lp{±iu) result from eqs. 
(4.101) and (4.133), we can be sure that these do not have singularities in the upper 
half plane, JRcu > 0. Also we assume, that these functions take finite values in a; = 
and for a; — >■ oo. 

The integration in (B.l) is converging for cu ^ and for a; — >■ oo due to the 
cancelations which take place in the parentheses. We would like to consider both 
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contributions separately. For the corresponding integrals to converge we introduce 
regularizations on both ends and represent 



i lim 

r->0 
R—i-oo 



/ — (p{iuj) - / — (p{-iuj) 

Jr ^ Jr ^ 



(B.2) 



In the first integral we change the integration variable, u —cu. Further we add and 
subtract integrals along pathes 7^ and 7^. Both are half circles in the upper half of 
the imaginary plane with radii r and R such that 



7 = [-R, -r] U 7, U [r, R] U -fn 
is a closed path. Accordingly, we rewrite the integrals in the form 

f doo , . , f dco . . . f dcu 
/ — (p[—tu!)— / — (p[—tuj)— / — (p[—tuj) 



i lim 

r-^0 



(B.3) 



(B.4) 



Now, for the functions we are considering, the integral over the closed path is zero. 
From the other two integrals, in the limit, we get contributions from the half poles. 



7 = 7r((^(0)-(^(oo)), 



(B.5) 



which is the desired formula. 
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